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ABSTRACT 

This paper presents ionization-corrected measurements of the carbon abundance in intergalactic 
gas at 4.0 < z < 4.5, using spectra of three bright quasars obtained with the MIKE spectrograph on 
Magellan. By measuring the C iv strength in a sample of 131 discrete H i-selected quasar absorbers 
with p/p > 1.6, we derive a median carbon abundance of [C/H]= —3.55, with lognormal scatter 
of approximately a ~ 0.8 dex. This median value is a factor of two to three lower than similar 
measurements made at z ~ 2.4 using C iv and O VI. The strength of evolution is modestly dependent 
on the choice of UV background spectrum used to make ionization corrections, although our detection 
of an abundance evolution is generally robust with respect to this model uncertainty. We present a 
framework for analyzing the effects of spatial fluctuations in the UV ionizing background at frequencies 
relevant for C iv production. We also explore the effects of reduced flux between 3-4 Rydbergs (as 
from He n Lyman series absorption) on our abundance estimates. At He n line absorption levels 
similar to published estimates the effects are very small, although a larger optical depth could reduce 
the strength of the abundance evolution. Our results imply that ~ 50% of the heavy elements seen 
in the IGM at z ~ 2.4 were deposited in the 1.3 Gyr between z ~ 4.3 and z ~ 2.4. The total implied 
mass flux of carbon into the Ly-a forest would constitute ~ 30% of the IMF-weighted carbon yield 
from known star forming populations over this period. 

Subject headings: quasars : absorption lines; intergalactic medium 



1. INTRODUCTION 

Intergalactic C iv has been studied in the spectra of 
quasars at essentially all redshifts, from the very local 
neigh borhood in the observed-frame UV (jCooksev et al.l 
[20101) to redsh ift six in the observed-frame near-infrared 
(iBecker et al.l [20091: iRyan- Weber et al.l [20091: iSimcod 
120061: iRvan- Weber et al.ll2006D . Yet most of our detailed 
knowledge about the ionization-adjusted metal content 
of intergalactic matter comes from st udies at z ~ 3 
(jSchave et al.l 120031: ISimcoe et al.ll2004D . In this interval 
the observed- frame transitions of both C iv and Ly-a fall 
in an optimal region for optical spectrometers, in terms 
of both sensitivity and lack of foreground contaminants. 

However, even at z ~ 3 the density of intergalactic 
gas is sufficiently low, and the UV ionizing background 
field is sufficiently intense, that most carbon is ionized to 
the C v state and higher. At typical IGM densities for 
z ~ 3, about 2% of all carbon atoms are in the C IV state. 
This fact, coupled with the low value of both the particle 
density and the overall heavy element abundance, leads 
to an extremely weak C iv signal in the more tenuous 
regions of the cosmic web. Consequently studies of metal 
abundances have required extremely high signal-to-noise 
ratio data (~ 50 - 200). 

To date the most sensitive surveys have either utilized 
statistical methods to tease signal fr om the distribution 
of C iv and H i pixel optical dep ths ( Sch aye et alJl2003l 
lEllison et al]l2000b lSongailal[200l . or used the O VI line 
which has a higher ionization potent ial (and therefore a 
stronger signal) to trace metallicity (jSimcoe etalll200l 

1 This paper includes data gathered with the 6.5 meter Mag- 
ellan Telescopes located at Las Campanas Observatory, Chile 

2 MIT-Kavli Center for Astrophysics and Space Research 

3 Alfred P. Sloan Research Fellow 



iDave et al.lll998HBergeron et al.ll2002fl . When the same 
UV backgrounds are used to estimate ionization correc- 
tions, both the C IV pixel method and the O VI method 
yield an IGM abundance of roughly [C/H]= —3.1 or 
[0/H]= -2.8 at z ~ 3. 

The O VI line does not lend well to studies of abun- 
dance evolution at high redshifts, because it is blended in 
the Ly-a forest, whose overall opacity increases at higher 
z. Above z ~ 3 it becomes effectively impossible to de- 
blend O vi doublets from H I forest lines, even in high 
resolution spectra. C iv is seen at these redshifts, with 
a distrib ution of streng ths that is nearly independent of 
r edshift (ISongai all2005h . 

iSchave et al.l ( 20031) systematically examined the ion- 
ization corrected carbon abundance between 3.0 < z < 
4.0, finding no statistically significant evidence of evolu- 
tion. However the bulk of this sample was centered at 
lower redshift, with only one object at > 4. 

In this paper, we revisit the question of abundance 
evolution by obtaining three new high-SNR spectra of 
z ~ 4.5 quasars. We use these to make a detailed com- 
parison of C iv absorption at z ~ 4.3 with the O vi and 
C iv measurements we obtained earlier at z ~ 2.5. 

The z ~ 4.0 — 4.5 range is a "sweet-spot" for observing 
C iv because it is the interval when the fraction of car- 
bon in the triply ionized state peaks over cosmic time. 
As one moves to higher redshifts, the predominant ioniza- 
tion state moves to progressively lower levels due to the 
combined (1 + z) 3 increase in the baryon density and de- 
cline of the hard-UV radiation field provided by quasars. 
Thus a survey of individual C iv systems at z ~ 4.3 is 
comparably sensitive to a survey for O vi at z ~ 2.5 or 
a statistical detection of C iv at z ~ 3. These different 
surveys all probe gas structures with baryonic overden- 
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Figure 1. MIKE spectra of the C IV region for the three objects used for the z ~ 4.3 sample. At top, the mean telluric absorption 
spectrum determined by standard star observations is shown for reference. The masked region at 8350Ain the spectrum of BR0353 is at 
the location of a known Mg II interloper. 

sities of p/p ~ 1.5 relative to the cosmic mean at their 
respective epochs. 

In Section [2] we describe the new observations and 
methods for data processing; Section[3]describes our sam- 
ple selection, and Sectio n [4] describes our observational 
measurements. Section 15.11 discusses ionization correc- 
tions applied to derive the abundance measurements pre- 
sented in Section [5] The resultant cosmic abundance 
evolution and implications for galaxy formation are dis- 
cussed in Section [7] Throughout we assume a cosmology 
with n M = 0.3, 57a = 0.7, and H = 71 km/s/Mpc. 
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2. OBSERVATIONS 

We observed the three bright, southern quasars listed 
in Table 1 using Magellan/MIKE over several different 
observing runs between 2004-2006. Spectra were ob- 
tained using a 0.875"slit, which yields a combined ve- 
locity resolution of 14km s _1 , as measured from concur- 
rent exposures of ThAr arc lamps. Extraction was per- 
formed using a custom ized software package for MIKE 
(|Bernstein et alj|2003f) . which makes use of the Poisson- 
limited 2D sky subtraction methods outlined in IKelsonl 
(|2003D . The Echelle orders of individual exposures were 
flux calibrated using observations of hot spectrophoto- 
metric standard stars. Finally, the calibrated orders 
from all exposures for each object were combined onto 
a single 1-D wavelength grid and converted to vacuum- 
heliocentric units. 

At z = 4.0 — 4.5, the C IV transition falls between 
A =7745A and 8520A. This region is rich with telluric 
absorption features from the atmospheric "A" and "B" 



bands of molecular oxygen and water. We carefully re- 
moved these features from our QSO spectra by fitting 
transmission functions to observations of our standard 
stars interspersed throughout each night. 

The final 1-D spectra are shown in Figure [TJ The re- 
sultant signal-to-noise ratio in the C IV region ranges 
from 20 to 40; in the Ly-a forest the ratio is somewhat 
lower but still suitable for measuring H I given the much 
greater strength of the features being measured. 

3. SAMPLE SELECTION 

To choose our spectral search regions for C iv, we 
began by fitting Voigt profiles to the entire H I Ly- 
a forest using the vpf it0 software package. We selected 



H i 



> 10 14 ' 5 as the minimum threshold column den- 
sity for an absorber to b e included in the C IV ca ndidate 
sample. For a standard lHaardt fc Madaul (|2001[ ) model 
of the UV background field, an absorber with this column 
density and [G/H]=-3.0 would have N c IV = Kpcm" 2 — 
roughly the detection limit of our data. 

4 http:/ /www. ast.cam.ac.uk/~rfc/vpfit. html 
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Figure 2. Montage of absorbers and corresponding best-fit Voigt profiles, illustrating the use of partial profiles for multiple Lyman series 
transitions. For each plot, the column density and redshift of the line shown in red are given at top. Blue curves show blended absorption 
for other lines with : > 10 14 ' 5 , which also are included in the C IV sample. 



At z ~ 4.3, this H I column density corresponds to a 
baryon ic overdensity of p/p — 1.6 relative to the cosmic 
mean (|Schavel [20011 ). Coincidentally, this threshold is 
nearly identical to the minimum p/p where individual 
O vi lines can be detected at z ~ 2.5; C iv lines can only 
be seen at higher overdensity at z ~ 2.5. Throughout the 
present paper we quote comparisons to both O vi and 
C iv samples at lower redshift. 

3.1. Reliability Tests for iVu ; Fits 

When fitting for JVh i at z > 4, much of the Lyman 
series is accessible provided there are no strong Lyman 
limit systems. We used at least three transitions of the 
Lyman series, and often 4 or 5 to constrain the column 
densities of each H I line. If no unsaturated Lyman series 
line or profile wing was available, a line was discarded 



from the sample. 

Figure [2] shows several examples of H I absorbers fit in 
this manner, selected to span the range of column den- 
sities in our sample. Clearly there is significant blending 
of lines in the Ly-a and higher order transitions, and it is 
rare to have a single absorber with many unblended tran- 
sitions. However roughly 60% of the sample has at least 
one Lyman series transition completely free of blending, 
and an additional 20% were significantly constrained by 
the unsaturated wings of blended lines in the high or- 
der transitions. Figure [3] shows a graphical summary of 
which transitions were used to constrain each line in the 
sample, over the full range in redshift. 

In practice, very few systems were rejected because of 
saturation in all Lyman series lines. For each sightline 
we truncated our C iv search below the redshift where 
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Figure 3. Illustration of which Lyman series lines were used in 
H I fits at varying redshifts. For each system, multiple transitions 
are used. Horizontal lines indicate components with at least par- 
tial saturation (which sometimes contain information in absorption 
wings). Open circles are transitions which may be blended but 
are not saturated in absorption and hence provide upper limits on 
.Nhi- Solid circles indicate distinct, apparently unblended lines 
which align in velocity with other Lyman series absorption. 

the Ly-7 line becomes saturated from Lyman limit (or 
completely blended Lyman series) absorption. The high 
redshift boundary for each sightline is located 5000 km 
s _1 blueward of the QSO's emission redshift. 

In total, our H 1 sample contains 139 lines, with 13 ad- 
ditional candidates discarded from the sample because of 
saturation in the Lyman series. Just one of the discarded 
systems is associated with a possible C iv absorber. 

When using high order Lyman transitions to constrain 
iVfji, there is danger of over-estimating the column be- 
cause of blended Ly-a absorption at lower redshift con- 
taminating the signal. To assess this potential error, we 
generated Monte Carlo realizations of the Ly-a forest in 
our redshift range to test whether the fitting procedure 
introduces systematic errors. 

The details of Ly-a forest line samples are not exten- 
sively characterized above z ~ 3.5 and it is beyond our 
present goals to provide a comprehensive statistical anal- 
ysis of the H 1 forest. However, as a byproduct of our fit- 
ting procedure, we obtain a determination of the H I col- 
umn density distribution at z ~ 4.2, shown in Figure |4l 
As at lower redshift, it follows a power-law form with 
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Figure 4. The H I Column Density Distribution Function, 
/(jv H I ,X) as measured in three sightlines at 2 ~ 4.3 (black points), 
and five sightlines at z ~ 2.4 (red points). The distributions are 
offset by roughly 0.5 dex, or a factor of three in normalization. 
Solid line shows a AT -1 power law, which was used as input for 
the Monte Carlo calculations described below. 
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Figure 5. Input vs. recovered column densities in our Monte 
Carlo tests of the vpf it method for selecting the H I sample. We 
detect an overall scatter of ~ 0.1 dex in determining iVjj 1, with 
systematic offset of 0.04 dex. Ionization corrections render the 
effective error in [C/H] smaller than the error in Nn 1 by a factor 
of 3. 

incompleteness roll off at Ah r < 10 14 cm~ 2 — wel l below 
the cutoff A H , of our C iv sample. iKim et all (|2001f > 
analyzed Ly-a forest spectra at z = 2.31 and z = 3.35, 
finding /(Ah T ) oc N^\ A6 . This same shape provides 
an accurate fit to our observed /(A) if a slightly higher 
overall normalization is used to account for the higher 
density of forest lines at z ~ 4.2. 

For the Monte Carlo tests, we assumed this form 
for /(A), with an overall normalization evolving as 
(1 + z) 217 ([Kim et alj 12001) to approximate the evolu- 
tion in density of lines. Artificial spectra were generated 
using lines drawn from these distributions, with the full 
Lyman series included for each absorber. Ly-a lines at 
lower redshift were intentionally added into the Ly-/3 and 
higher order regions at appropriate density to simulate 
blending and contamination. Finally, Gaussian noise was 
injected into each spectrum with an amplitude set by the 
noise arrays in our sample spectra. We do not claim the 
resulting spectra to be a perfect representation of the 
Ly-a forest at z > 4, but they form a useful tool for 
bootstrapping estimates of our systematic errors in Ah i 
measurement from blending. 

We ran VPFIT on our simulated spectra without 
knowledge of the input parameters, in identical fashion 
to the true data set. Figure [5] shows a comparison of the 
best-fit Ah i determinations with the true input values, 
as a function of Ah i and redshift. We find a scatter of 
±0.1 dex between the true and fitted values, with a small 
systematic (median) offset of 0.04 dex in the sense that 
the fitted Ah i over-estimate the true values. Individ- 
ual outliers can miss by 0.5-0.7 dex in either direction, 
but [65,85]% of systems are determined to better than 
[0.1,0.2] dex. 

3.2. Comparison Sample 

To study evolutionary trends, we compare the z ~ 4.3 
C iv sample with a separate sample of abs orbers at z ~ 
2.4, t aken from previously published work (|Simcoe et al.l 
I2004D . In this work, O VI column densities were measured 
for a sample of 230 absorbers with Am > 10 136 cm~ 2 . 
A smaller sample of 83 systems were used to measure 
C iv absorption at A H 1 > 10 14 cm -2 . 

In Figure [6l we plot the cumulative distribution of p/ p 
for each sample. We estimate the overdensity at each 
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Figure 6. Distribution of overdensity in our H I -selected samples 
at z ~ 4.3 and z ~ 2.4. The C IV sample at 2 ~ 4.3 probes smaller 
overdensities than the z ~ 2.4 C IV sample, by a factor of 2. The 
distribution is very similar to O VI samples at lower redshift. 

H I ab sorber location using the formalism of lSchave et al.l 

UH - 10 Cm U-3 x 10" cm-* ) Ti ^{oIg) 

Here, T4 represents the temperature in units of 10 4 K, Ti2 
represents the flux of ionizing photons at the H 1 ioniza- 
tion edge, and f g represents the gas to total mass frac- 
tion in Ly-ct forest systems, which is assumed to be near 
fifc/fijw. To obtain overdensity, we normalize this den- 
sity by the mean gas density, fibp c (l + z) 3 . We assume 
T4 to be constant a nd unity over the z = 2.4 — 4.3 range 
(|Schave et al.|[2000h . 

Figure[5]shows that the C IV sample at z = 4.3 probes a 
median overdensity that is 2 x lower than the C iv sample 
at z ~ 2.4, despite having a limiting TVh 1 that is 3x 
higher. The O vi sample at lower redshift probes almost 
an identical range in gas overdensity to the high redshift 
C iv sample. Therefore, comparing C iv at both redshifts 
removes any ambiguities in the [C/O] relative abundance 
but means that the high redshift sample is probing lower 
densities. On the other hand, comparing [O/H] at z ~ 
2.4 and [C/H] at z ~ 4.3 traces the abundance evolution 
at fixed overdensity for an assumed value of [O/C]. 

Unfortunately this quantity is degenerate with one's 
choice of UV backgrou nd spectrum used for calculating 
ionization corrections (ISimcoe et alJ I2004t iSchave et all 
12001 lAguirre et alj 12008ft . In the discussions below we 
assume [0/C]> in the IGM, consistent with ioniza- 
tion by a soft spectrum compos ed of both galaxies and 
quasars (jHaardt fc Madau1 l2001). This model is favored 
by pixel opt ical depth measurements of C Hl /C IV and 
O vi /Si iv (jSchave et al.ll2003l : lAguirre et al.ll2008ft . 

We caution however that observations of the most 
metal-poor damped Ly-a systems show a recovery to- 
ward [O/C]=0 at very low metallicities (|Pettini et al.l 
12008ft . Solar [O/C] ratios can be achieved with a hard 
background composed of quasar light with negligible 
galaxy contribution. 

3.3. C iv Sample Definition and Measurements 

For each of the H 1 absorbers in the sample, we ex- 
amined the corresponding wavelength of the C IV dou- 
blet. Using VPFIT, we fit C iv absorption profiles to any 
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Figure 7. Cumulative distribution of C IV column densities at 
z ~ 4.3 and z ~ 2.4. 

absorbers within a velocity range of ±20km s _1 of the 
H I centroid. In blended systems, we matched individual 
C iv components to their corresponding H I components 
where possible. If a single H I line contained multiple 
components of C iv then the total column density for 
these components was recorded. Finally, if no C iv is 
detected we used VPFIT to calculate ler upper limits on 
the C iv column density, using the error array. 

The final sample therefore consists of a mix of de- 
tected systems with measured column densities, and non- 
detections with upper limits. We used a > 3cr cut as the 
threshold between detections and non-detections. For 
the non-detections, VPFIT usually finds a formal solu- 
tion Nc iv ± crc iv with Ac iv < 3trc iv! the upper limit 
is then Nc i V + 3oc iv 

The final sample ofz = 4.0-4.5Hi and C iv redshifts 
and column densities is presented in Table 2. The three 
sightlines contain 131 systems with A^h i > 10 14 5 Among 
this sample, 39 lines have C iv detections, while 92 are 
non-detections. The approximate limiting C iv column 
density is between 1 — 2 x 10 12 cm~ 2 . 

4. CUMULATIVE DISTRIBUTIONS OF C iv , AND 

N c iv/JVh i 

Figure [7] shows the distribution of C iv column den- 
sities in samples drawn at the two different redshifts 
considered. Because the C iv samples contain a mix- 
ture of detections and upper limits, we have plotted the 
cumulative distribution as determined by the Kaplan- 
Meier estimator. The Kaplan- Meier product limit is 
a fundamental tool of survival statistics, which deals 
with censored datasets su ch as the one presented here 
(|Feigelson fc Nelson! 11985ft . Its application for measur- 
ing abundan c es fro m absorption lines is discussed in 
ISimcoe et aD (|2004ft . 

For the Kaplan-Meier distribution to accurately repre- 
sent the true distribution, two conditions must be sat- 
isfied. First, the upper limits must be statistically in- 
dependent of one another, which is clearly true for the 
resolved, discrete absorption lines in our sample. 

Second, the probability of a non-detection must not 
depend on the criterion used to select the sample. Since 
we are measuring C iv detections and upper limits for a 
sample selected on the basis of H I column density, this 
condition should hold. There may be slight discrepancies 
at the very high end of the H I sample because of local 
chemical enrichment. But because the H I column den- 
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Figure 8. Cumulative distribution of the C IV /H I ratio at z ~ 4.3 
and z ~ 2.4. This is the fundamental observed quantity in our 
analysis. 

sity distribution is a steep power law, the vast majority 
of our sample is at the low end where the censoring only 
depends on variation in signal-to-noise in the C iv region 
and is therefore quite random. 

Returning to Figure it is clear that the distribution 
of C iv evolves very little between z ~ 2.4 and z ~ 4.3. 
The fact that the cumulative distributions only reach ~ 
40% stems from the survival analysis — recall that C iv is 
only detected in 39 of 131 systems, or 30%. Most of the 
upper limits are clustered near low Nq iv , but some are 
at slightly higher values. Kaplan-Meier accounts for the 
difference between the fraction of detections (30%) and 
the value of the cumulative distribution near the limiting 
N Clv (40%). 

The constancy of the C IV distribution wit h redshift 
has be en noted by other authors, most notably ISongailal 
( 2005) . Our result is consistent with these previous find- 
ings. 

For studies of the IGM metallicity, our fundamental 
observable is Nq iv/^Vh i- Because C iv evolves very little 
with redshift but the H I opacity increases, we expect 
the C iv /H I ratio to decrease toward higher redshift. 
The cumulative distribution of this quantity is shown 
in Figure [H with the expected difference between z ~ 
2.4 and 4.3. Note that this is the distribution of the 
Nc iv/^Yhi ratio for each system in the sample, not the 
ratio of the C iv and H I distributions. The shapes are 
very similar at the two epochs, with the difference coming 
from a systematic shift downward by roughly 0.4 dex, or 
a factor of 2.5. 

5. IONIZATION CORRECTIONS 

To translate the evolution in Nq iv/^h i hito an evo- 
lution in the carbon abundance, we apply an ionization 
correction to each individual absorber as follows: 




Here, fc iv an d /h i represent the triply ionized frac- 
tion of carbon, and neutral fraction of hydrogen, re- 
spectively. The last term is a zero-point offset to nor- 
malize relative to the solar abundance. Throughout 
the paper we have assumed t he m eteoric Solar abun- 
dances of iGrevesse &: Sauvall (|1998[) . where A car bon = 
8.52 on a scale with Androgen = 12. Several revisions 
have been proposed to the Grevesse solar abundances 



Figure 9. Mean spectrum of the UV background radiation 
field, at z ~ 4.3 and z ~ 2.4. We show the calculations of 
both Haardt & Madau (2001, including a galaxy contribution) and 
IFaucher-Giguere et al.l (120091 ) for compariso n. The spectra are very 
simila r at high redshift; at lower redshift the Fauchcr-Gigucrc et al. 
(2009) version is harder, leading to higher overall abundances and 
a correspondingly larger evolution in the observed [C/H] ratio. 

(jAllende Prieto et al.l [20021 12001 IBolwegerl 12001 ; we 
have kept them for consistency with prior work by many 
authors on IGM abundances. 

Ionization fractions for carb on and hydrogen we re de- 
termined using the CLOUDY ([Ferland et alJ[l99l soft- 
ware package. We assumed the IGM is optically thin 
(consistent with the low H 1 column densities in the sam- 
ple), and in photoionization equilibrium. 

5.1. Comments on choice of the UV background 
spectrum 

The largest uncertainty in our ionization correction 
arises from the choice of a particular prescription for the 
ionizing background spectrum. 

Several groups have recently provided improved con- 
straints on the overall intensity of the background field 
at 1 Rydberg - denoted as Jgi2- This can be derived from 
the opacity of the Ly-a forest, which relates in turn to 
the H 1 photoionization rate Tw In a sam ple of 63 high 
resolu tion spectra spanning 2.0 < z < 5.5, Becke r et alJ 
(2007) find T 12 = [0.8,0.5] at z = [2.4,4.3], where 
by convention Ti2 is n ormali zed to units of 10~ 12 s _1 . 
IFaucher-Giguere et alJ (|2008af > find T 12 = [0.6, 0.4] over 
the same range using slightly different analysis tech- 
niques on a set of 86 hi g h res olutio n spectra. Earlier 
analyses by iBolton et all (|2005l ) and iScott et all (|2000[ ) 
derived slightly higher estimates closer to Ti2 = 1. How- 
ever, all of these studies find that the evolution of the 
H I photoionization rate is quite weak over the redshift 
range studied here, perhaps 10 — 30% smaller at z = 4.3 
but not much more. Throughout the rest of the paper, 
we use the values from iBecker et all (|200l . For a power 
law ionizing spectrum oc i/~ 1,8 , this corresponds to a flux 
at the H 1 ionization edge of log(Jgi2) = [—21.5, —21.7] 
at z = [2.4,4.3]. 

Even if the normalization of the background spectrum 
at 1 Rydberg does not evolve, it is generally agreed that 
the shape of the spectrum will change. The quasar lu- 
minosity function fa lls off between z — 4 and z = 2 
(Hop kins et aHl2007L While ionizing photons at z = 2.4 
come from both quasars and galaxies, the background at 
z = 4.3 contains a higher fractional co ntribution from 
galaxies (jFaucher-Giguere et alJ l2008af h This should 
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Figure 10. Neutral fraction of hydrogen as calculated by 
CLOUDY, compared with the analytic approximation of Equation 
3. This relation is redshift-independent, except for the implicit de- 
pendence on the IGM temperature and Ti2. It has a logarithmic 
slope of ^ dex change in /h : per unit dex in Nn : . 

lead to a softer UV background at earlier epochs, with 
fewer photons at the higher energies which regulate the 
ionization balance of heavy elements like carbon and oxy- 
gen. 

We have tested two main prescriptions of the UV back- 
ground with our data: t he commonly used version from 
lHaardt fe Madaul (120011). and an indep e ndent new de- 
termination from IFaucher-Giguere etaLl rf2009h . In Fig- 
ure [9] we show these spectra, with the energies of rel- 
evant transitions for the C iv ionization balance la- 
beled. The Haardt & Madau spectrum includes con- 
tributions from both QSOs a nd galaxies; the QS O por- 
tion is determined from the iCroom et al.l (|2004l) lumi- 
nosity function, while the galaxies are added through 
population synthesis modeling, assuming a 10% escape 
fraction of photons beyond the Lyman limit. The 
Haardtfc Madau spectra have been normalized to match 
the observational constraint s on F12 described above. 
IFaucher-Giguere et al.l (I2009D t ake t he QSO luminosity 
function from iHopkins et al.l ((2007) which falls more 
steeply toward higher redshift, and add a galaxy con- 
tribution to keep the evolution of Ti2 constant. 

When the Haardt & Madau spectrum is renormalized 
in CLOUDY to match observations of Ti2 at 13.6 eV, 
its soft X-ray background is also artificially boosted by a 
factor of ~ 10 at energies above 10 Rydbergs. This has 
consequences for the C iv ionization balance because it 
affects the C v to C vi transition. 

In their calculations, Haardt & Madau treat the UV 
and X-ray backgrounds independently, trying to match 
each with observed luminosity functions by balancing the 
relative importance of Type I and II QSOs. In its raw 
form, the Haardt & Madau X-ray background is slightly 
lower at z=4.3 than at 2.4 (as would be expected), but 
after renormalization to matched the observed T12 at 
1 Ry, the situation reverses, such that the X-ray back- 
ground is higher at z = 4.3 than it is at z — 2.4. Since 
the X-rays originate in AGN which are proportionately 
fewer at high redshift, this situation is probably unphys- 
ical. So, we applied a downward correction of 0.8 dex 
to the z — 4.3 HM spectrum above 10 Ryd to soften 
it by a comparable amount as the UV background and 
bring it back into agreement with the original X-ray in- 
tensity. Coincidentally, this downward correction brings 
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Figur e 11. CLOUDY carbon ionization balance at redshift 4.3 
for the Haardt & Madau ( 2001) model of the UV background (top) 
and Fauchcr-Gigucrc et al. (2009, bottom) models. The histogram 
at top indicates the cumulative distribution of lines in our sample, 
with the minimum, median, and 90th percentile shown by vertical 
lines. At this redshift, the C IV fraction peaks in the IGM, at 
fc iv ~ 0.5. 

the HM spectrum into fairly close agreement with the 
IFaucher-Giguere et al.l (|2009t ) spectrum, which by con- 
struction does not suffer from this effect of renormaliza- 
tion. 

While the HM and IFaucher-Giguere etaLl ((2009) spec- 
tra do not agree in every detail, their general shapes 
match fairly well, particularly at z = 4.3. At lower red- 
shift the difference is larger owing to assumptions on the 
AGN contribution as well as the treatment of He II reion- 
ization. However, in both models the main change from 
z = 2.4 to z = 4.3 is a factor of 10-20 decrease in the 
hard UV background. It should be noted that neither 
m odel captures the He II "sawtooth" effect described 
in iMadau fe Haardt (2009); this is discussed in Section 

MM 

5.2. H I ionization fractions 

Since the Ly-a forest is optically thin and in pho- 
toionization equilibrium, a calculation of /h i is relatively 
straightforward and may be made even without detailed 
knowledge of the spectral shape. In equilibrium, the ion- 
ization and recombination rates may be balanced as: 

/h in H T « n 2 H R (3) 

where the recombination rate R ~ 4 x lO -13 ^ 0,76 
cm s _1 , uh is given in Equation 1, and we assume the 
gas is mostly ionized. Equations 1 and 3 may be com- 
bined to derive the hydrogen ionization fraction as a func- 
tion of (observed) H I column density, plotted in Fig- 
ure [10] alongside its numerical CLOUDY determination. 
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Figure 12. CLOUDY carbon ionization balance at z ~ 2.4, again 
for both models of the UV background. The mean H I column 
density is lower, as is the ionization fraction in C IV. The lower 
ionization fraction may be directly mapped to the inferred change 
in abundance. 

This relation is independent of redshift, except implic- 
itly through the (weak) evolution of the ionization rate 
r, and the temperature T4. 

5.3. C IV ionization fractions 

The ionization balance of C iv is a far more difficult 
calculation, because it is governed by photons at frequen- 
cies where the UV background spectrum is less well con- 
strained observationally. Figure fTTI shows the CLOUDY- 
derived ionization balance for carbon as a function of 
H 1 column density, at z = 4.3 (the balance at z = 2.4 
is shown in Figure [T2l . To calculate the ionization frac- 
tions, we first estimated the total hydrogen density asso- 
ciated with each H 1 column density as calculated from 
Equation 1. This value of nn was input into CLOUDY 
along with the prescriptions shown in Figure |9] for the 
shape and amplitude of the background spectrum. To- 
gether, the determinations of the density and ionizing 
background flux determine the ionization parameter U 
for the cloud which allows CLOUDY to calculate the 
abundance fraction of each carbon ionization state. 

The top panel of each figure shows the cumulative dis- 
tribution of H 1 column densities for the C iv samples, as 
a visual aid to indicate what ionization corrections are 
appropriate for which fractions of the sample. Vertical 
lines on the bottom panels indicate the sample minimum 
Vh i, and also its median value and 90th percentile. Plots 
are shown both for the Haardt & Madau and Fauchcr- 
Giguere forms of the background spectrum. 

These plots show that the C iv ionization state pre- 
dominates at z ~ 4.3, at a value of fc iv = 0.50 — 0.65 
for the median H I system in the sample. This range re- 



flect s the difference between thelHaardt fc Madaul (|2001[ ) 
and IFaucher-Giguere et al.1 (|2009f) backgrounds, respec- 
tively. Roughly equal amounts of carbon are in higher 
and lower ionization states. 

Conversely, the median system at z ~ 2.4 contains only 
4-10% C iv , with the vast majority of carbon inhabiting 
higher ionization states. So, the typical ionization cor- 
rection at z — 4.3 is an upwards factor of 1.5-2.0 (0.2-0.3 
dex), compared to a correction factor of 10 — 25 (1.0-1.4 
dex) at lower redshift. 

6. ABUNDANCE RESULTS 

6.1. Evolution from z = 4.3 to z = 2.4: Qualitative 
Results 

Before presenting the full carbon abundance distri- 
bution results, it is instructive to examine the evolu- 
tion of the median system using qualitative arguments. 
From Figure [5J the median H I absorber decreases in its 
C iv:H I ratio by 0.4 dex between z — 4.3 and z — 2.4. 

Suppose now that there were no evolution in the [C /H] 
ratio. To achieve this, our measured 0.4 dex decrease in 
C iv /H I would need to be offset by a corresponding 0.4 
dex increase in log (/h i//c iv), according to Equation 2. 

We have seen already that /h i depends on Nn , but 
only very weakly on redshift. For a median log(-/Vn i) = 
14.5 at z = 2.4 and 15.2 at z = 4.3, the H I ion- 
ization fractions can be read from Figure I10| yielding 
log(/H 1,4.3) ~ log(/H 1,2.4) ~ 0.5. 

Thus for fixed fc iv, the first term in Equation 2 is 0.4 
dex lower at z ~ 4.3, while the second is 0.5 dex higher. 
The first term is observationally determined, and the sec- 
ond is fairly model-independent, and they approximately 
cancel each other out. The degree of evolution in [C/H] 
therefore corresponds almost exactly to the change in 
fc iv with redshift. 

Even if we take the most conservative possible hypoth- 
esis - that the UV background spectrum at z = 4.3 is 
identical to the spectrum at z — 2.4 - then U will de- 
crease because of the increase in gas density at higher 
redshift. But observational studies of the QSO and 
galaxy luminosity functions at these redshifts indicate 
that the background spectrum should soften as it is in- 
creasingly dominated by galax ies toward higher redshift 
(jFaucher-Giguere et al.ll20u8al) . These factors cause an 
increase in fc iv which leads to a decrease in implied 
[C/H]. Thus the observed distribution of Nav/Nm, 
coupled with minimal assumptions about ionization, sug- 
gest that [C/H] is lower at z = 4.3 than it is at z = 2.4. 

To estimate by how mu ch, consider the CLOUDY cal- 
culations from both the lHaardt fc Madaul ()200lD and 
IFaucher-Giguere et al.l ()2009f ) spectra. These illustrate 
that the C IV ionization fraction for the median system 
increases by 0.7 dex toward higher redshift. Taken to- 
gether, these qualitative arguments suggest that [C/H] 
in the median system is lower at z ~ 4.3 by 0.5-0.6 dex 
(a factor of ~ 3) than it is at z = 2.4. 

They also suggest that the scatter in the C iv to H 1 ra- 
tio may be an interesting diagnostic of fluctuations in 
the hard UV background. The scatter in this ratio re- 
flects the convolved scatter in the intrinsic [C/H] ratio, 
measurement errors, and variance in the ionization cor- 
rection. Even if one makes no assumptions about the 
relative weighting of these three factors, one may set up- 
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Figure 13. Kaplan-Meier cumulative probability function for 
carbon abundance in the IGM at z ~ 4.3 and z ~ 2.4. Note 
that the low redshift curves sample the IGM at a median over- 
density roughly 2x larger than the high redshift sample. Solid 
curves indicate lognormal probability distributions. At z ~ 4.3, 
the curve represents a median abundance of [C/H]=-3.65, and 
a = .8. At lower redshift the median rises to [C/H]=-3.1 
for the IHaardt fc MacTaul 120011) background, or [C/H]=-2.8 for 
IFauche r-Giguere ct al. (2009), with similar scatter. The difference 
of 0.5-0.7 dex represents a 3 — 5x increase in the median abundance 
for the low redshift sample. 

per bounds on the scatter in the hard background under 
the conservative assumption of zero measurement error 
and abundance scatter. 

6.2. The Cumulative [C/H] Distribution at z ~ 4.3 

Figure Q2] shows the Kaplan-Meier distribution of 
[C/H] at z = 4.3 and z = 2.4, derived by applying ion- 
ization corrections individually to each C iv system (or 
its upper limit). The distri butio ns are shown for both 
the IHaardt fc Madaul (|2001[) and IFaucher-Giguere et al.l 
(2009]) for ms of the UV background radiation field. Re- 
call that we show the [C/H] distributions for each red- 
shift to eliminate ambiguity in [O/C], but the z = 4.3 
sample probes gas with a median overdensity of p/p « 3, 
whereas the C iv sample at z = 2.4 only reaches a me- 
dian p/p w 6 (the sample minima are p/p w 1.6 and 
p/p~3). 

At z — 4.3, our sample reaches nearly to the median 
[C/H] (50% on the plot), which lies near [C/H]= -3.55, 
or ~ 1/3500 the Solar abundance. Roughly half of the 
lines in our sample have 3tr upper limits restricting them 
below this value, so the shape of the distribution below 
the median is not known. Above the median, the distri- 
bution appears to be very crudely lognormal, with the 
upper quartile falling above [C/H]= —3.1. The distri- 
butions are very similar for the two choices of UV back- 
ground radiation spectrum at z = 4.3. The smooth solid 
curve shows the cumulative distribution for a lognormal 
probability distribution. Our best fit lognormal model 
shows a median abundance of [C/H] =-3.55, and lognor- 
mal deviation of a = 0.8 dex. 

At z = 2.4, we find a medi an abundance of [C/H]= 
-3.1 with a = 0.8 dex for the IHaardt fc Madaul (pOOlf ) 
form of th e background s pectru m, consistent with prior 
results of ISimcoe et al.l (120041) and close to, t hough 
slightly higher than those of ISchave et al.l (|2003[h The 
IFaucher-Giguere et al.l (|2009f ) form of the background 
spectrum results in a slightly higher carbon abundance 
but similar scatter, at median [C/H] =—2.8 and a = 0.7 
dex. The higher median results from the increased flux in 



the IFaucher-Giguere et al.l (|2009l ) spectrum near 5 Ryd- 
bergs, at the C IV to C v ionization edge. This flux 
enhancement favors higher ionization and hence larger 
ionization corrections. 

Comparing the two distributions, we see that the 
sample median carbon abundance has increased by ~ 
0.45 — 0.65 dex (depending on the background model) 
during the epoch between z — 4.3 and z — 2.4, al- 
though the density probed at z ~ 2.4 is a factor of 2 
larger. It is interesting to compare this result to that of 
ISchave et al.l d2003h . who find d/dz([C/R]) = +0.08±g;?§ 

and d/d(log(/9))([C/H]) = 0.65^ ^. Not accounting for 
density, our measurement implies a derivative of roughly 
—0.26 dex per unit redshift, a factor of 3.3 higher than re- 
ported by Schaye et al. Moreover the evolutionary trend 
has the opposite sign as Schaye's; our new result evolves 
in the expected direction (i.e. increasing abundance with 
time), although the prior result was basically consistent 
with zero evolution. 

Since Schaye et al. reported weak correlation with red- 
shift but strong correlation with density, it is worth con- 
sidering whether our measured evolution is simply an 
artifact of the different densities probed at the different 
redshifts. According to the regression analysis of Schaye 
et al., a 2x change in density alone could account for 
« 0.3 dex difference in [C/H] between the samples, but 
their reported evolution of +0.08 dex per unit redshift for 
Az = 1.9 would counteract the density effect by 0.15 dex, 
yielding a total change in [C/H] of 0.15 dex to our 0.5. 
Schaye et al interpret the small statistical significance of 
their evolution coefficient and its non-intuitive sign as a 
non-detection. We have detected a larger difference in 
[C/H] between our high and low redshift samples, in the 
expected direction, even accounting for the difference in 
density. 

To decouple the effects of density and true evolution, 
we ran our analysis for a subset of our high redshift 
sample restricted to log(iVH i) > 10 14 7 , thus ensuring 
that [C /H] is measured at fixed overdensity between the 
two redshifts. The resulting Kaplan-Meier distribution 
is shifted to higher [C/H] by 0.2 dex for this factor of 2 
increase in p/p, suggesting that the abundance is indeed 
lower in regions of smaller overdensity, but by a smaller 
degree than measured by Schaye et al. If we budget 
0.2 dex of the total 0.5 change seen in Figure Q2] to the 
change in gas density, the remaining 0.3 dex represents 
the temporal evolution signal. 

As described earlier, we may also compare our mea- 
surements of [C/H] at z ~ 4.3 to pre-existing measure- 
ments of [O/H] at z ~ 2.4, which cover an identical range 
of overdensity bust must be corrected for any non-Solar 
value of [O/C]. We estimate the value of [O/C] by com- 
paring the [O/H] and [C/H] distributions at z ~ 2.4, and 
applying this same value at z = 4.3 (i.e. assuming no 
evolution in relative solar abundances). Figure fT4l shows 
the result of this calculation. The [C/H] distribution is 
shown in the thick solid line, while the [O/H] distribution 
is the thin line; here the [O/H] curve has been shifted 
downward by 0.25 dex, simulating the pr edicted [C/H] 
for a best-fit [O/C]=0.25. Other studies (|Simcoe et al.l 
l200i ISchave et al.112003 : lAguirre et af1 l2008) have gener- 
ally found positive [O/C] distributions for soft photoion- 
izing spectra, with values near [0/C]~ 0.5. Our slightly 
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Figure 14. Kaplan-Meier Distribution function for Carbon and 
Oxygen in the z ~ 2.4 IGM, used for calibrating the [O/C] ratio, 
which is best fit with [O/C]=0.25. A positive [O/C] has been found 
by many authors, though the value used here is slightly smaller 
than other estimates in the literature (see text). 
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Figure 15. Kaplan-Meier Distribution function for [C/H] at z ~ 
4.3 and [O/H]-0.25 at z ~ 2.4. Subject to estimates of [O/C], this 
distribution is our best representation of the change in [C/H ] at 
fixed overdensity. The curve at z ~ 4.3 is the same as in Figure Tl3l 
the lognormal fits shown here for z ~ 2.4 have a median [C/H]=- 
3.1, -2.9 for the two models of the UV background. 

smaller value likely results from a slightly higher normal- 
ization of t he background spectrum to match rece nt ob- 
servations ([Becker et ail I2007t iBolton et all 12005ft. The 
best-fit [O/C] for the (iFaucher-Gieuere et al J 12009ft UV 
background is near [O/C]=0, which is to be expected 
considering that this sp ectrum is much harder than the 
(|Haardt fc Madadl200lft background at z = 2.4. 

Figure [15] shows the distributions of [C/H] at 
z = 4.3, [O/H1-0 .25 at z = 2.4 for the 
IHaardt & Madaul (l200lft U V background, and [O/H] for 
the lFaucher-Giguere et ail ([2009ft model (which has best- 
fit [O/C]=0). Lognormal curves are show for each model. 
The two z = 2.4 models have median abundances of —2.9 
and —3.1, with similar scatter. Evidently there is a 0.5 
dex evolution in the median [C/H] at fixed overdensity 
for [O/C]=0.25; adopting [O/C]=0.5 would reduce the 
absolute change by 0.25 dex. 

Figure [16] shows this evolution as evidenced in the indi- 
vidual sample measurements at each redshift. The large 
circles show our median estima tes at each redshift, while 
the central point is taken from lSchave et al.1 (|2003ft 

Taken all together, these calculations indicate that the 
intergalactic carbon abundance increased by 0.3-0.5 dex 
between z ~ 4.3 and z ~ 2.4, with little change in the log- 
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Figure 16. Illustration of the individual measurements for each 
line in the full sample. Large red dots indicate the median abun- 
dance at each redshift, calculated using the Kaplan-Meier product 
limit estimator. The median points lie below the center of the 
scatter because of the inclusion of numerous upper limits. The 
solid p oint at z = 3.1 represents the measurement of lSchave et al.l 

normal abundance scatter. The range quoted reflects un- 
certainty in various choices for the UV background, the 
relation between density and metallicity, and the value 
of [O/C] used. 

A change of 0.3-0.5 dex represents a factor of ~ 2 — 3 
increase in the median intergalactic carbon abundance 
at fixed overdensity. It suggests that roughly half of the 
carbon seen in the IGM at z ~ 2.4 was deposited in the 
1.3 Gyr interval between z ~ 4.3 and z ~ 2.4. 

6.3. Uncertainties 
6.3.1. Ionizing Background Fluctuations 

The results presented above all assume a spatially uni- 
form shape and intensity of the ionizing background, fol- 
lowing the convention of all similar previous studies. Al- 
though we find similar results for two representations of 
the mean background (Haardt & Madau, and Faucher- 
Giguere), we did not explore the uncertainty or scatter 
that would be introduced by place-to-place variations in 
the radiation field. 

This subject has largely been ignored in the previous 
literature on IGM abundances, but recent work on the 
rcionization of H I and He II contains arguments that are 
relevant to observations of heav y element lin es. 

Following the develop ment of[Fardal ct al. (1998j) and 
iMeiksin fc Wh ite (200$, the variance in the background 
intensity at a given energy is a function of the source 
density of photon emitters at that energy, the mean free 
path of photons with that wavelength through the IGM, 
and the variance in source spectrum from object to ob- 
ject. When ionizing sources are rare, differ in SED from 
object to object, and photons have a small mean free 
path, the radiation seen by a random absorber tends to 
be dominated by a single source and will vary signifi- 
cantly from place to place. 

The production of H I ionizing photons at z = 4.3 
occurs mostly in galaxies which are relatively numerous, 
and the mean free path for these photons is large because 
H I is already reionized. So, many sources will be con- 
tained within an H I photo-attenuation volume and the 
background should be comparatively uniform at 1 Ryd 
(except in the immediate vicinity of QSOs). 

Overall variations in [C/H] are therefore more sen- 
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sitive to the C iv ionization balance, whose atomic 
transitions are resonant with photons between 3 — 
30 Rydbergs. These photons originate in less nu- 
merous AGN, which could lead to an increased vari- 
ation in spatial intensity. Since our high redshift 
sample is measure d at the probable l eading edge of 
He II reionization ([Reimers et al.lll997|: [j akobsc n et al. 
19941 iHeap et al.ll2000l; iDayidsen et alj|1996t iKriss et al. 
200lt iFechner et all 120061: iShull et al.l 120101: iSongaila 



1998D . the IGM opacity from photoelectric absorption 
must be considered above the He n ionization potential 
(4Ryd). 

For photons of en ergy E, the mean free path is given 
as (jFurlanettoll2009l ) 



A = 0.41 



5.3 



E 



4 Ryd 



He 



onMpC (4) 



At z = 4.3, 

ii ~ 1 since He II reionization is only in 
its earliest phases. For C iv — >C v ionizing photons with 
4.74 Ryd, this implies A < 1 Mpc because of He n at- 
tenuation, while for C v — >C vi (24.2 Ryd) the path is 
90 Mpc. Thus the most important effect we must con- 
sider is stochastic fluctuation in the 4.74 Ryd background 
governing the C iv to C v transition. 

At z = 4.3, simulations indicate that 10-20% 
of the universe resides in ionized "bubbles" of 
He in (|McQuinn et al.l 120091) . He n Ly-a observations 
focus only on the interiors of these bubbles, where the 
ionized fraction is high and He n absorption troughs are 
not totally saturated. The lower oscillator strength of 
C iv renders it observable in both He ill ionized bub- 
bles, and the He n -neutral walls in between. In the 
bubble walls, He II absorption strongly attenuates the 
UV background near 4 Ryd, affecting the ionization bal- 
ance of carbon. Inside of bubbles, stochastic variations in 
the proximity of hard photon sources and attenuation by 
He II Lyman-limit systems generates intensity variations 
both above and below the mean spectrum. In the follow- 
ing two sections we consider each of these environments 
in separate detail. 

6.3.2. UV background fluctuations in the interior of 
He in bubbles 

To study ionization variations of carb on inside 
He III bubbles, we follow the arguments of iFurlanettol 
(2009), who develops a simple Monte Carlo framework 
suited to the problem. The method was originally in- 
tended for studying variations in the 4 Ryd background 
in He ill bubbles at z ~ 3; a minor set of modifications 
renders it useful for studying the C iv — >C IV transition 
at at E — 4.74 Ryd and z = 4.3. Essentially all opaci- 
ties in the model are scaled down to reflect the smaller 
cross section of He II at higher energy, and fluxes from 
the ionizing sources is also scaled down to reflect their 
SED. 

We defer a full description of the method to Furlan- 
etto's paper. Very briefly, all quasars are assumed to live 
inside of He ni-ionized bubbles. The space outside of 
bubbles is opaque to He Ii-edge photons; inside of bub- 
bles there is still attenuation from the He n analog of Ly- 
man limit systems, but between these systems the mean 
free path is large. 

Because QSOs are rare, the typical bubble contains 
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Figure 17. Distribution function of flux at the C IV — >C v ion- 
ization edge, in He Ill-ionized bubbles at z ~ 4.3. The flux is 
normalized in units of the mean, (J). Scatter above the mean is 
dominated by proximity to individual QSOs, which enhances the 
background over a wide range of frequencies above 1 Ryd. Scatter 
on the low side arises from stochastic variations in the density of 
nearby QSOs and shielding from He II patches. 

only one to a few sources, so a Monte Carlo approach 
is used. Ionizing sour ces are drawn at ran dom from the 
luminosity function of I Hopkins et al.1 (|2007t ). and placed 
uniformly within each bubble volume (no sources are lo- 
cated outside of bubbles) . Then, the spatial variation of 
the background J is calculated by compiling trial statis- 
tics with varying QSO luminosities, spatial locations in 
the bubble, Lyman limit attenuation, and UV spectral 
indices from the sources. 

The Monte Carlo calculation outputs a distribution of 
j = J I (J), where (J) represents the mean value of the 
background in a fully He il-ionized IGM. The J distri- 
bution at the C iv — ;>C V ionization edge can be used to 
examine variations in our C IV ionization corrections. 

The spread of j depends strongly on the radii of the 
He III bubbles, denoted as R, since variations diminish 
once bubbles be come large enough for abso rbers to "see" 
multip le QSOs. iMcQuinn et al.l (|2009H and lBolton et al.l 
( 2009) have shown in simulations that the distribution of 
R does not evolve strongly with redshift and is governed 
primarily by QSO duty cycle. Typical values range from 
20 — 40 comoving Mpc, with the spread at any redshift 
larger than the difference between z = 4 and z = 2. The 
spread in j also depends on the photoattenuation length 
ro, which parameterizes the distance that C iv ionizing 
photons travel from their p oint of origin befo re reaching 
a He II lyman limit system. IFurlanettol (|2009| ) calculates 
ro for He n ionizing photons at z = 2,3,4 and demon- 
strates that its luminosity-weighted mean evolves very 
little. C iv ionizing photons have slightly higher energies; 
their attenuation is still dominated by He II absorption, 
but the relevant ro increases very slightly because of the 
smaller He II cross section at 4.74 Ryd. 

Figure [171 shows the distribution of J at 4.74 Ryd in 
logarithmic (top panel) and linear (bottom panel) forms, 
for bubbles ranging from R = 15 to 50 comovi ng Mpc. In 
all cas es we have assumed r = 35 Mpc, as in Furlanettol 
(2009). As expected, the shape is very similar to what 
was derived for He n bv lFurlanettol (|2009l ). except the dis- 
tributions are slightly narrower because of the increased 
transparency of He n at 4.74 Ryd. Above the mean 
background flux, the curves converge for all bubble sizes. 
IFurlanettol (|2009f) observed this same effect for He II; it 
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Figure 18. Model UV background spectra used to evaluate the ef- 
fects of variations on our inferred fC/H] estimates. T he thick solid 
line represents the mean [Haardt & Madau| 11200 II ) background. 
Curves below this mimic the suppression of 4.74 Ryd photons by 
He II Lyman limit systems (inside of bubbles) or the neutral walls 
between bubbles. Curves above the mean represent the proximity 
zones of QSOs. 

represents environments where the background is dom- 
inated by a single source. In this case the variance is 
driven by the shape of the QSO luminosity function and 
the random distribution of QSO-absorber spacings. 

According to Figure [TZl variations on the high-side are 
limited to a factor of ~ 10 x enhancement in the back- 
ground. Because this enhancement is driven by proxim- 
ity to a single QSO, it should lead to an increase in the 
UV background across a wide range in frequency. When 
examining the C iv ionization balance in these neigh- 
borhoods, we generate an "enhanced" toy model of the 
background spectrum which boosts the flux for all ener- 
gies above 1 Ryd by a constant factor of 5 — 10. Below 
1 Ryd the background remains the same, since it should 
still be dominated by integrated starlight from more nu- 
merous galaxies. 

Below the mean J, the distribution in Figure[17]is more 
extended and varies with bubble size. The increased vari- 
ance for small bubble sizes results primarily from shot 
noise in the number of quasars: by R = 15 Mpc most 
bubbles have either no source or one source, and sources 
outside the bubble are shielded from view. This effect 
is strongest at the He n ionization edge, and still signifi- 
cant at 4.74 Ryd. But the bubble walls are optically thin 
to photons with E < 4 Ryd and also for higher energy 
photons including the 24.2 Ryd C v — >C vi edge. Thus 
it is appropriate simply to use the mean UV background 
spectrum for H I , as well as all carbon transitions ex- 
cept C iv — v (and possibly C in — >C iv , see Section 
I6±5l) . 

To model effects at 4.74 Ryd, we attenuate the 
lHaardt fc~ Madau (200l|) spectrum by an absorbing col- 
umn with r = To(-E/4Ryd)~ 3 . This produces the desired 
suppression at 4.74 Ryd while preserving the cosmic av- 
eraged spectrum in optically-thin regions. The normal- 
ization t is adjusted to produce a factor of 10-100 de- 
crease in the background at 4.74 Ryd, consistent with 
the distribution in Figure 1171 

6.3.3. UV background fluctuations in the He n -neutral 
inter-bubble medium 

In simulations of He II reionization, He in bub- 
bles only permeate 10-20% of the IGM by volume, so 



Figure 19. CLOUDY calculations showing correction factors to 
the ionization term in Equation 2. In regions where the 4.74 Ryd 
background is suppressed, little or no change is apparent. In re- 
gions of enhanced flux and high density, we may overestimate the 
abundance by up to ~ 1 dex. 

we must consider the shape of the spectrum in the 
He ii walls between bubb l es. T his scenario is discussed 
in iFaucher-Giguere et al.l (|2009fl : essentially flux at the 
He II edge is completely suppressed because in He n- 
neutral regions the optical depth can easily exceed r ~ 
100. Once again, the background at lower energies is spa- 
tially uniform, being optically thin after H I reionization 
(and, for low energies, dominated by galaxies). Likewise 
at higher energies the He n cross section declines and the 
hard flux recovers to a uniform background. 

We therefore model the background spectrum in 
He II walls the same as in the "suppres s ed flu x" regions 
of ionized bubbles: a lHaardt fc Madaul (|2001D standard 
mean background having He n absorption blueward of 
the edge with r = ro(-E/4Ryd) -3 . However in the neu- 
tral walls the normalization of optical depth is much 
higher. Whereas attenuation from Lyman limit systems 
in bubbles reduces the flux by a factor of 10-100, the 
completely neutral bubble walls with r ^ 100 experience 
attenuation by 44 orders of magnitude. 

6.3.4. The effect of spatial fluctuations on [C/H] estimates 

Figure [18] shows the range of ionizing background 
models we considered when investigating spatial fluc- 
tuations. The thick solid line represents the mean 
lHaardt fc Madaul ([2001) quasar+galaxy spectrum at z = 
4.3. Four models with He n attenuation of to = 
1, 3, 10, 100 fall below the mean. Smaller values of To rep- 
resent conditions on the lower half of the J distribution 
inside of He in bubbles, the larger values are characteris- 
tic of the He n bubble walls. Above the mean curve are 
two models with enhanced hard UV flux (i.e. proximate 
to QSOs), shown in dot-dashed lines. 

We ran each of these spectra through our CLOUDY 
model grid to examine the change to the ionization frac- 
tion as a function of column density. Figure [T5] shows the 
result of this exercise, where we plot the change in ion- 
ization correction for each different choice of background 
spectrum. A value of zero indicates no change in [C/H] 
for a given line; a positive value implies a lower metal- 
licity for a given set of C iv and H I measurements. 

The solid lines near zero represent the four "attenu- 
ated" models with varying optical depth at 4 Ryd, in- 
creasing upwards. For these models the H I ionization 
fraction does not change at all because the universe is al- 



13 



_5 jp- — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — ■— 
: Effect of Hell Lyman Series Sawtooth 

? q 4 l MH2009 sawtooth on v~ us Background 

X> '- 

: As above + all Lyman series r doubled 

§ 0.3: 



0.0 L 
1 4.0 



15.0 
log(N HI ) 



Figure 20. Effect of flux suppression at the C III ionization edge 
from the He II Ly-ct forest "sawtooth" spectral modulation. 

ready optically thin to 912 A photons at z ~ 4.3 and the 
background is dominated by flux from numerous galax- 
ies. The C iv ionization fraction does not change either, 
which results in a very small total correction of < 0.2 
dex to [C/H] across our range of Nh r . 

The small change in fc rv is surprising at first glance 
given that the flux at the C iv — »C v ionization edge 
changes by many orders of magnitude. The reason this 
effect is so small is that even for the mean spectrum, 
over 50% of the carbon is already in the C iv state. So, 
a reduction in the C iv — >C V rate — which would increase 
fc iv — cannot change the C iv fraction by more than a 
factor of ~ 2. 

A much larger correction is obtained when considering 
upward fluctuations in the UV background from proxim- 
ity to bright quasars. While this change can be signif- 
icant for selected systems, the number of affected data 
points should be small. Major upward fluctuations occur 
only in the largest bubbles surrounding rare, luminous 
QSOs, and even in these bubbles the probability of ly- 
ing in a proximate region is low according to Figure 1171 
Such an error would push our abundances lower, enhanc- 
ing the evolutionary trend discussed in Sections 16.11 and 
IO 



6.3.5. He II Sawtooth 

iMadau fc Haardtl (|2009l ) have recently emphasized the 
importance of a "sawtooth" imprint left on the ionizing 
background spectrum between 3-4 Ryd, from line ab- 
sorption of the He n "Lyman" series. This relates to the 
C iv balance because the transition from C in to C iv oc- 
curs at 3.5 Ryd. These authors' updated model of the 
UV background including the sawtooth modulation leads 
to a decrease in flux at the C in ionization edge of 0.2- 
0.3 dex at z — 3. They also explored models where an 
artificially delayed He II reionization led to a reduction 
in flux at the C in edge of over a full dex. This latter 
case may be more appropriate for studying the z = 4.3 
IGM since He II reionization should not be complete at 
this epoch and the resulting opacity may be quite high. 
Statistical studies of metal li nes at high redshift pro vide 
some evidence for this effect (|Agafonova et al.l 120051 ) 

For the mean spectrum with no He II sawtooth, the ma- 
jority of carbon is in the C IV state at z ~ 4.3. Therefore, 
in principle a softening of the background should shift the 
balance toward C in and C II, reducing fc iv arL d increas- 
ing [C/H] estimates in the process (Equation 2), thereby 
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[C/H] (Abundance) 

Figure 21. C omparison of Kap l an-Me ier distributions for the 
standard z = 4.3 Haardt & Madau (2001) model, the z = 2.4 stan- 
dard model, and the z = 4.3 model incorporating a 1 dex decre- 
ment at 3-4 Ryd from the sawtooth modulation. This modification 
changes the shape of the Kaplan-Meier distribution, such that the 
high [C/H] regions do not evolve, but the low [C/H] regions do. 

reducing the evolution signal. 

To explore this effect, we obtained a copy of the saw- 
tooth spectrum from F. Haardt, and used it to recalcu- 
late ionization corrections. In its raw form, the sawtooth 
model has a very hard spectrum that is disfavored by 
IGM abundance studies at lower redshift ()Schave et al.l 
l2003tlSimcoe et aT1l2004[ ). It produces an abundance gra- 
dient that decreases with density, leaving very high abun- 
dances in weak Ly-a forest lines. 

Rather than use this form as-is, we instead examined 
the ratio of the sawtooth spectrum (HM+S in their Fig- 
ure 1) to a separate model calculated with identical in- 
put parameters but no sawtooth included (HM in their 
Figure 1). This isolates the effect of the He II Lyman se- 
ries, which we then multiply into the softer background 
spectra shown in Figure GO Clearly this is not a self- 
consistent way to model He n absorption in the IGM, and 
IMadau fc Haardtl (|2009f) point out that even their models 
do not fully capture the patchy nature of He II reioniza- 
tion that could be important here. But lacking a full sim- 
ulation suite, our approach captures the spectral shape 
of the sawtooth, and provides an heuristic tool for exam- 
ining its impact on our results. 

Figure [201 illustrates the change in [C/H] resulting from 
use of our modified sawtooth spectrum, as a function of 
AThi- The solid line shows the MH09 model; for the 
dashed curve we arbitrarily increase the optical depth 
across the sawtooth by a factor of 2. Depending on the 
degree by which the spectrum is changed, the correction 
can be anywhere from 0.08 to 0.2 dex, and always in the 
upward sense. The sawtooth modulation therefore will 
decrease the evolution signal between our z = 4.3 and 
z — 2.4 samples. 

Figure ED shows the Kaplan-Meier distribution in 
[C/H] obtained by running our survival analysis with 
the HM+S background at z = 4.3. The two compar- 
ison curves show the distributions at z = 2.4 and 4.3 
for the unmodulated spectrum. For our modified HM+S 
spectrum, the [C/H] distribution shifts uniformly to the 
right by ~ 0.06 dex; a similar plot using the background 
with 2x inflated He II optical depths would shift ^0.15 
dex higher than the distribution from the unmodulated 
spectrum. 

Recall that we observe a 0.5 dex change in [C/H] be- 
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Figure 22. Si IV /C IV ratio for systems in the sample with de- 
tected C IV . Model curves show the predicted trends for vari- 
ous degrees of sawtooth absorption in the UV background spec- 
trum. The solid c u rve in corporates no sawtooth, dotted line is the 
IMadau fc Haardt (200S|) form modified as described in the text. 
The dashed and dot-dash lines are for sawtooth spectra with opti- 
cal depth arbitrarily multiplied by 2 X and 4 X , respectively. These 
models, which would be needed to substantially shift the [C/H] 
distribution, are disfavored by the few measurements shown here. 

tween z — 4.3 and z = 2.4; if one forces [O/C]=0.5, or 
requires a strong correlation between density and metal- 
licity, the evolution could be as small as 0.3 dex. A very 
strong sawtooth modulation in the UV background could 
in principle remove 0.15 dex of our evolution signal. If 
one applies both the sawtooth spectrum and forces abun- 
dances to [O/C]=0.5, it is possible to derive a solution 
with only 0.10-0.15 dex of evolution in the carbon abun- 
dance, much closer to the no-evolution scenario. This 
situation arises only for a set of assumptions stressed to 
a particular end of our parameter space, but its possi- 
blility should not be ignored. 

The chief effect of the 3-4 Ryd decrement at z ~ 4.3 is 
to take carbon from the C IV state and move it to C III, 
so we should in principle be able to see an increase in 
the C ni/C iv ratio using the C in 977A line, or even an 
increase m C n/C iv at 1334A. We searched our data for 
systems where this measurement could be made, but un- 
fortunately the results do not provide strong constraints 
on the magnitude of the sawtooth effect. The C in lines 
tend to be blended with Ly-a forest absorption, and 
C II lines are too weak to detect in all but our strongest 
few systems. One of the strongest constraints may come 
from pixel-optical-depth analyses. The C ni/C iv ratio is 
not a good d i agnos tic at z > 4. However, the models in 
ISchave eT al. (2003) — which show no trend of [C/H] with 
rcdshift — may yield a declining cosmic abundance from 
high values in the past toward low values at present us- 
ing a sawtooth background — a solution we are generally 
prejudiced against. 

The Si iv /C iv ratio is another diagnostic of spec- 
tral hardness which is accessible for a limited num- 
ber of systems in our sample. Figure [35] shows the 
calculated values, or upper limits (where no Si iv is 
detected). The solid lines show predicted values of 
Asi w/Nc iv for different prescriptions of the sawtooth, 
calculat ed using CLOUDY with Solar relative abun- 
dances (jGrevesse &; Sauval I1998D . At high Ah i where 
the most reliable measurements are located, the data are 
broadly consis tent with the non -s awtoo th spectrum and 
our modified Mad au fc Haardti (|2009l ) sawtooth form. 
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Figure 23. Error introduced in ffm, and hence [C/H], caused 
by errors in our estimate of the QSO continuum in the Ly-a forest. 
No systematic trend is apparent, and most of the scatter is at the 
~ ±1 dex level. 

Artificially boosted sawtooth absorption overpredicts the 
amount of Si iv seen in these systems, although the num- 
ber of points is too small to draw broad conclusions. This 
diagnostic suffers a slight degree of ambiguity because of 
possible non-solar values of [Si/C]. However many ab- 
sorbers at high redshift show [Si/C]> 0, in which case 
there would be an even stronger conflict for the heavily 
absorbed sawtooth spectra. 

6.3.6. Continuum Fitting 

At high redshifts an additional possible bias is intro- 
duced from uncertainties in the absolute continuum flux 
from the background quasar over the Ly-a forest region. 
We have followed the typical procedure of estimating con- 
tinuum levels by fitting a low-order spline across rela- 
tively unobscured regions of the forest. This procedure 
becomes increasingly inaccurate at high redshifts, as the 
opacity of the Ly-a forest increases, and it becomes in- 
creasingly difficult to find regions with transmission of 
u nity. 

IFaucher-Giguere et al.l (|2008bD present a treatment of 
this same problem in the context of Ly-a forest opacity 
measurements. Drawing mock spectra from Ly-a forest 
simulations at various redshifts, they performed manual 
continuum fits similar to the ones described here and 
examined the resultant deviation from the known, true 
continuum as a function of redshift. Their results indi- 
cate that at z — 4.2 (our sample redshift) a manual esti- 
mate could systematically place the continuum too low, 
by ~ 17%. By comparison, at z = 2.5, where O VI and 
other studies focused, the correction is of order 2%. 

To examine this effect, we created a "continuum 
corrected" version of our highest quality sample spec- 
trum (BR0353-3820) and refit H I column densities 
to the ful l Ly-q forest. Follow i ng the empirical ap- 
proach of IFaucher-Giguere et al.l ((2008b), we adjusted 
the continuum-normalized flux at each wavelength in the 
Ly-a forest downward by a factor of (1 — C{z)), where 
C(z) — 1.58 x 10~ 5 (1 + z) 5 - 63 , to account approximately 
for the bias introduced by incorrect continuum fits. We 
then re-calculated Ah i at each redshift for our full sam- 
ple. 

Figure |2"31 shows the results of the re-fitting exercise for 
BR0353-3820, where as before we use multiple Lyman 
series transitions to constrain Am- Somewhat surpris- 
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ingly, we found that the column density for the median 
system was nearly unchanged, or at least was offset by an 
amount smaller than the ~ ±0.2 dex point-to-point scat- 
ter. This may result from the fact that in each system, 
the JVh i measurement is driven by that Lyman series 
line which is closest to saturation but not bottomed-out 
across the profile. The fractional error introduced by 
poor continuum fits is smallest at low flux levels so by 
choosing the right transitions this is minimized. At least 
it seems clear that changes in the continuum level at 
z ~ 4 lead to errors at the level of < 0.5 dex in Ah i- 

In general, continuum errors enter into [C/H] in two 
ways. One is a simple error in Nq iv/Ah i. However, be- 
cause fc iv oc N^/f (Figure [ID|). we also underestimate 
the neutral fraction, which partly counteracts the mises- 
timate of Nc rv/Njj i. Thus the error in [C/H] is only ^ 
the logarithmic error in Nn i — in our case, ~ 0.15 dex of 
random scatter with no systematic offset. 

There is an additional, much smaller contribution from 
an incorrect estimate of fc ivj which is derived from 
Nn 1 . Over the column density range of our sample, the 
C IV fraction is near its peak. This results in a broad min- 
imum in the logarit hmic derivative of the C IV f raction 
with A H j. For the IFaucher- Giguere et al.l (|2009t l back- 
ground (as a representative example; see Figure [Til . 
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from 14.5 < log(A ? n i) < 16.5, with the derivative near 
zero between 15.0 and 15.75. The C iv column density 
in Equation 2 is not sensitive to continuum fitting errors 
and remains constant. 

Taken all together, the corrections to Ahi,/hi, and 
fc iv would lead to a combined correction in [C/H] of 
only 0.032-0.034 dex for each 0.1 dex change in A^hi- 
Moreover the sense of the change is such that the mea- 
sured abundance at z = 4.3 becomes lower as Nh i is 
increased, which tends to enhance the trend of declining 
carbon abundance with increasing redshift. 

6.4. Summary of systematic uncertainties 

To synthesize the results of this section: the main 
systematic sources of uncertainty in our abundance es- 
timates come from the UV background used to cal- 
culate ionization corrections, and the determination of 
H I column densities in the thick Ly-a forest at z ~ 
4.3. We have tried two choic es for the ionizin g back- 
ground: the commonly-used lHaardt fc Madaul (|200lD 
model with quasars + galaxie s, and a newer calculation of 
IFaucher-Giguere et al.l (|2009f ). Both versions of the spec- 
trum yie ld an evolving carbon abunda nce with redshift, 
with the IFaucher-Giguere et al.l (2009) version giving a 
slightly larger derivative. 

We investigated the effect of spatial variations in the 
UV background due to photoelectric absorption from 
He II , and stochastic fluctuations in the ionizing source 
population. These variations do not affect the H I neu- 
tral fraction, or transitions relevant to C iv other than 
the C iv to C v transition at 4.74 Ryd. Near this edge, 
the background can vary by factors of 10-100 inside of 
He ill-ionized bubbles, and many orders of magnitude in 
the neutral bubble walls. However, CLOUDY tests re- 
vealed that changes in the 4.74 Ryd flux, absent changes 



at other wavelengths, affects our abundances only at the 
0.1 dex level. It is possible that we have overestimated 
abundances for systems in the proximity zones of bright 
quasars by ~ 1 dex or more, but these systems are quite 
rare and should not drastically alter our distributions. 

The inclusion of "sawtooth" absorption from the 
He II Lyman series at 3-4 Ryd leads to a systematic in- 
crease of 0.06 dex in our median [C/H] for an absorption 
streng th consistent with the results of iMadau fc HaardU 
(2009) . If the sawtooth absorption strength is artificially 
inflated this offset becomes correspondingly larger. 

Abundance errors from a systematic underestimate of 
Nn i (as would be found from continuum fitting errors) 
are limited to < 0.1 dex, because for every 0.1 dex in- 
crease in Ah i, there is a corresponding increase of 0.07 
dex in /h i offsetting the effect, and little change in the 
C iv ionization fraction. 

Taken together these tests suggest that [C/H] abun- 
dances at z ~ 4.3 are less sensitive to systematic errors 
than one might naively expect. Different choices of the 
mean background will change the overall median abun- 
dance but do not affect our conclusion that the carbon 
abundance is evolving except under stressed sets of as- 
sumptions. Most of the effects studied above tend to en- 
hance the evolutionary signal, if we have underestimated 
their importance. 

7. DISCUSSION 

7.1. The Integrated Carbon Production Between 
z = 4.25 and z = 2.4 

The the total mass flux of carbon into the Ly-a forest 
may be estimated simply from our data by integrating 
the abundance- weighed H I column d ensity distribution. 
As described in iSimcoe et al.l {2004), the contribution 
of carbon atoms (in all ionization states) to the closure 
density may be estimated as: 



n c = l(| 
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where p c is the (current) closure density, and (C/H)q is 
the solar carbon abundance by number. The H I column 
density distribution /(jv H i,x) is defined as d 2 Af /dXdN-R : 
with Af being the number of absorbers in a survey of 
comoving pathlengh dX (Shown in Figure |4|) . 

The mean carbon abundance, shown in angled brack- 
ets, may be derived from the distributions in Figure 
1131 Assuming a lognormal parameterization of the abun- 
dance, the mean may be determined as 
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For our measured median abundance of [C/H] = —3.55 
and scatter of a — 0.7 dex, this amounts to a mean car- 
bon abundance of 10~ 281 at z ~ 4.3. The same equation 
applied to z ~ 2.4 yields a mean abundance of 10~ 236 
at lower redshift. 

We substitute these values into Equation 8 to derive 
Q c = 2.7 x 10~ 8 at z = 4.3, and 4.8 x 10~ 8 at z ~ 2.4. 
This represents a factor of 1.7 increase over an interval 
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of just 1.3 Gyr. Put another way, roughly half of the car- 
bon observed in the Ly-a forest at z ~ 2.4 was deposited 
into the IGM after z ~ 4.3. Even in the conservative 
limit where 0.2 dex of the difference between our red- 
shift points is attributed to density effects, the fraction of 
"recently deposited" carbon would be about one-third of 
the total at z ~ 2.4. The rise in intergalactic abundance 
suggests a connection with the concurrent rise in the cos- 
mic star formation rate and/or the luminosity density of 
AGN in this same redshift interval. Many of these metals 
may reside in enriched regions near high redshift galaxies, 
given their relatively recent production. 

Outflows are of course commonly seen i n the spectra 
of star forming galaxies at these redsh i fts (|Pettini et al.l 
[20011: ISteidel et all 120101: IBielbv et all 12010ft and there 
is evidence of spatial association between strong 
C iv systems and UV-selected galaxies at z ~ 3 
(jAdelberger et all 120031: ISteidel et all 120101 ). Our mea- 
surement would seem to connect these correlations in the 
local galactic environment to the global growth of metal 
abundance in the IGM, in a temporal sense. 

Estimates of the stellar population ages for z ~ 3 
gala xies typically range b etween 100 — 200 Myr to ~ 1 
Gyr (Shaplcv etalll2005h . This means that the archety- 
pal lyman Break galaxy at z = 2.4 could be built from 
scratch within the redshift interval of our C iv survey. 
The winds seen in 1BG spectra travel at several hun- 
dred km/s, so even if these galaxies began driving to- 
tally unbound winds at birth, the wind material could 
only have traversed several hundred kpc and not deep 
into the IGM. 

The C iv column densities studied here (Nc iv ~ 
10 12 ) generally fall below the range where correla- 
tions with galax ies ar e strongest (Ng TV > 10 13 5 — 
Adelberger et all 120031: ISteidel et ail 120101 ) . Moreover 
Songailal ( 20061 ) studied similar C IV systems to the 
ones presented here, concluding that the weaker lines 
dominating absorption selected samples are too quies- 
cent to trace actively evolving wind bubbles. However 
they do trace densities that would be found in the fil- 
aments seen in cosmic simulations, where galaxies also 
reside. These filaments may be populated with moder- 
ately recent, cooled wind relics that would manifest as 
the weaker C iv systems seen in our survey and other 
sensitive C iv searches in the ly-a forest, although fur- 
ther theoretical work would be required to assess if this 
scenario holds up in detail. 

7.2. Comparisons with the Global Star Formation Rate 

For any moment in time, the difference between the 
intergalactic carbon growth rate and the yield-weighted 
star formation rate measures the global efficiency of feed- 
back processes at ejecting metals from nascent galax- 
ies. Observations of individual star forming galaxies 
in the local universe indicate that galactic winds can 
be significantly loaded in mass and especially in met- 
als (jMartin et all 12002). During periods of intense star 
formation, it is possible that all metals produced in a 
burst will escape; in more quiescent periods these same 
elements would settle back into the galaxy's ISM. At 
z = 2.5 — 4.5 individual galaxy measurements are dif- 
ficult, but using our C iv sample we may estimate a 
globally averaged efficiency of metal ejection. 



First, we calculate the total stellar mass formed over 
our redshift interval b y integrating the co smic star for- 
mation rate density of iReddv et all (|2008f h Over the in- 
terval 2.4 < z < 4.3 (At = 1.34 Gyr), this yields a total 
formed stellar mass of 1.6 x 1O 8 M0/Mpc 3 , or an aver- 
age of about 0.16 Mq Mpc~ 3 yr _1 . We then estimate the 
carbon output re sulting from this star form ation activity 
using yields from lWooslev fc Weaver] (|1995fFI . The IMF- 
weighted carbon yield, calculated assuming a Salpeter 
IMF from 0.5 to 40 Mq and for solar metallicity, is 
1.2 x 10~ 4 . The yield for progenitors with M < 10M© 
(which mostly still lie on the main sequence) is set to 
zero. Folding this into the star formation rate, we find 
that the total production of carbon from core collapse 
supernovae between z — 4.3 and z = 2.4 is 

AM C « 2.1 x 10 4 M o Mpc~ 3 (Produced in galaxies). 

(8) 

In the previous section, we found that the total growth 
of intergalactic carbon came to Aflg = 2.2 x 10~ 8 , or 

^c-i.gxio-^ 1 - (9) 

Multiplying by the critical density, we obtain the inte- 
grated volumetric increase in carbon mass: 

AM C = 7.1 x 1O 3 M Mpc" 3 (Added to IGM). (10) 

A simple comparison of Equations 8 and 10 shows that 
the rate of carbon enrichment in the IGM is roughly a 
factor of ~ 3 lower than the rate of carbon production in 
stars at this epoch. This analysis is subject to substan- 
tial uncertainties in calculating the SFRD, IGM carbon 
abundance, and yields. It is somewhat remarkable that 
in spite of this, the carbon production and pollution rates 
are of similar magnitude. Taken at face value, the result 
suggests that galaxies at z ~ 2 — 4 keep roughly 70% of 
the heavy elements they produce, and donate the other 
30% back to the IGM. 

There is precedent for this behavior at both low and 
high r edshift. At the high redshift end, iSimcoe et all 
(2004) estimated galaxy yields based on a closed-box 
chemical enrichment model of the ly-a forest. This anal- 
ysis relates to the present one, except that it is an in- 
tegral calculation whereas the present version is differ- 
ential. However the results are similar: the closed-box 
model requires that galaxies lose at least ~ 10% of their 
heavy element mass to enrich the IGM to observed levels 
by z ~ 2.4, and possibly more. 

Our escape fraction of ~ 30% is also comparable to 
those derived from very limited samples of dwarf star- 
bursts with outflows in the local universe. iMartin et all 
(2002) report ejection efficiencies of 15 — 20% for the 
dwarf starb urst NGC1569. For the l ocal p rototype star- 
burst M82, iStrickland fc Heckmanl (|2009l ) report a star 
formation rate of 4-6M Q yr" 1 with mass outflow rates of 
1-2 M Q yr" 1 . If the SFR is weighted by yield and the 
ejecta have abundances above yjj-Z© (observations sug- 

5 The IMF correction from UV luminosity to stellar mass is 
counterbalanced by the IMF-weighting of the yield, and it is the 
UV-bright stars which matter most for enrichment. As long as the 
same IMF is applied consistently for both calculations the detailed 
form matters less for metal production. We thank the referee for 
raising this point. 
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gest Z ~ 5Z Q in the wind), the metal outflow rate even 
exceeds the metal production rate, presumably because 
of mass loading from the ISM. 

While selected individual galaxies in the local universe 
have mass and metal outflow rates similar to these val- 
ues, the high redshift growth in abundance would require 
that essentially every galaxy drive outflows with this ef- 
ficiency. 

7.3. Metal Mixing and the Smoothing Scale for 
Abundance Measurements 

The abund a nces m easu red in this wo r k — as well as in 
iSimcoe et all (|200l and iSchave et all (|200l — are un- 
derstood to be smoothed on scales comparable to the 
Jeans length for the structures in question, which fall in 
the range ~ 50 — 150 kpc. Ly-a forest clouds have small 
overdensities relative to the cosmic mean, and are gen- 
erally thought t o arise in str uctures that are marginally 
Jeans unstable (|Schavell2001[ ). This sets the relation be- 
tween A^h i and nn seen in simulations, which we have 
exploited in calculating our ionization corrections. It in- 
dicates that the sizes, densities, and H I ionization frac- 
tions employed here should accurately represent the true 
physical conditions inside of H I absorbers. 

However the mixing of C iv on ~ 100 kpc scales is nei- 
ther well-mapped observationally or well-resolved in nu- 
merical simulations. The information that we do know 
suggests that C iv is not smoothly distributed through- 
out Ly-a forest clouds. This is suggested observation- 
ally by the slight velocity offsets often observed between 
H I and C iv lines at similar redshift. Moreover obser- 
vations of starburst outflows in the local universe show 
a highly non-uniform geometry. Typically there is a uni- 
form hot X-ray halo powered by the feedback source, 
but cooler gas that would be seen in absorption is often 
clumpy or filamentary. 

Several recent hydrodynamic codes have explicitly 
tracked the fate of metals flowing out of galaxies, and 
these also tend to find a very non-uniform distribution, 
with metal absorption coming from small, higher den- 
sity regions embedded in the more diffuse H I m edium 
(jCen fc Chisara2010tlQppenheimer fc Davel2006D . How- 
ever these simulations employ SPH or other particle 
based methods, tracking heavy element mixing by tag- 
ging outflow particles. These methods will underestimate 
the diffusion of metals on scales at or below that of a par- 
ticle, which are crucial to the mixing process. 

Several authors have studied the sizes of C IV absorbing 
regions along the lin e of sight using ionization modeling 
(|Simcoe et al.l l2004). or using transverse separa tion pro- 
vided by gravitational lenses ([Rauch et al.ll200ll) . Invari- 
ably these studies find small absorber sizes of ~ 1 kpc or 
less, sometimes much smaller. The e xtreme argument fo r 
non- uniform mixing is presented by (jSchave et al.|[2007| ). 
based on fitting heavy element components on the wings 
of larger C iv profiles. 

Using simple scalings, Schavc ct al. (2007) examine the 
fate of non-uniform, high metallicity clumps outside of 
galaxies. If these systems start out as hot clumps from 
superwinds, they are overpressurized with respect to the 
intergalactic surroundings, and will expand until they 
reach pressure and temperature equilibrium with the am- 
bient medium. Once in equilibrium, the metal-rich patch 
also reaches a common density with the IGM. This pro- 
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Figure 24. Comparison of the C IV column denisty distribution 
from sample spectra with th e redshift-independent power law de- 
termined by (Songaila 2001). For this plot only, we use Songaila's 
cosmology for direct comparison. Small number statistics set in at 
Nc iv ^ 10 14 , while incompleteness is apparent at JVc iv i; 10 12 5 . 
The few strong systems missed have minimal impact on the sur- 
vival analysis; the statistics are dominated by the more numerous, 
low column density absorbers. 

cess occurs on timescales of 1-10 Myr, essentially instan- 
taneously with respect to the Hubble time or even the 
star formation age of galaxies at z ~ 2 — 4. 

Small patches of high metallicity embedded in the Ly- 
a forest clouds may therefore dominate the C iv absorp- 
tion profile, but the H I profile is dominated by the ex- 
tended, metal poor Ly-a forest absorber. Because their 
densities are the same once equilibrium is reached, the 
ionization correction for C iv (which depends on nn de- 
rived from TVjji) should actually still be accurate for 
calculating fciv- This implies that we have correctly 
measured the relative amounts of hydrogen and carbon 
for calculating [C/H]. However this estimate is essen- 
tially averaged over the whole Ly-a absorber. The actual 
abundance in any single volume element may be either 
substantially lower or higher depending on the degree of 
mixing. 

7.4. Comparison with Previous Studies 

A number o f authors have studied the evolution of 
C iv at z > 2. Songaila] f|2005T ) notably evaluated flc iv 
and the C iv column density distribution function from 
z = 2.1 to z ~ 6. A key finding of this work was the lack 
of evolution observed in the C iv CDDF. Our results 
are broadly consistent with this result. Figure [7] illus- 
trates that our data yield a very similar column density 
distribution of C iv at redshifts 2.4 and 4.3. The Kaplan- 
Meier distribution is somewhat different than a standard 
CDDF in that we have not scaled by redshift pathlength, 
but by reporting as a percentage of the sample the dif- 
ferent pathlcngths of the samples normalize out. Also, 
the KM distribution better captures information about 
non-detections, which explains why C iv is detected in 
only about 50% of our sample (which in turn consists of 
stronger Ly-a forest lines). 

To test whether the C IV systems found in our sight- 
lines are consistent with other surveys from the litera- 
ture, we show the classic C iv CDDF in Figure 124"! made 
from a count of (C iv -sele cted) absorbers in our spec- 
tra. The power-law fit from iSongailal ()200lD is shown as 
a solid line, this fit agrees with our data over the range 
12.4 < log(ATc iv ) < 14.0. Both we and Songaila are sub- 
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stantially complete in the range 13 < log(Ac iv) < 14. 
A straight sum of systems in this Nc iv interval yields 
iv = 1-3 x 10 -8 for our sample at z = 4.3, while her 
result (scaled to SIm = 0.3) yields 1.1 x 10~ 8 , in good 
agreement. 

Using her ful l sample with 13.0 < log(Ac IV ) < 15.0, 
iSongailal ()2005h estimates that fl c iv = 1.8 - 2.7 x 1CT 8 
over our redshift range. In Section l7Tl we calculated the 
total ionization corrected Qc for an H I -selected sample, 
with result 2.2 x 1(T 8 and 4.8 x 10~ 8 at z = 4.3, 2.4, 
respectively. From Figure [TTJ we see that the C iv frac- 
tion in our sample ranges from 35 — 50% at z = 4.3 
to 17 — 36% at z = 2.4. These imply an approximate 
iv ~ 1.0 — 1.5 x 10~ 8 , on the low side of Songaila's 
range but generally consistent. It is not surprising that 
our estimate made in this alternate way is slightly low, 
since about half the signal in raw measurements of iv 
is produced by the rare, strong systems (i.e. outliers) 
that we have not included in our statistical calculation. 
The strong systems picked up in flc i V are not well sam- 
pled here and may reside in circum-galactic environments 
t hat are locally enrich ed. 

ISchave et all (|2003l ) performed a systematic investiga- 
tion of C iv at similar redshifts and, importantly, in- 
cluded ionization corrections in their analysis. Unlike 
our method, which relies on measurements of line col- 
umn densities and upper limits. ISchave et a l. (2003) per- 
formed a pixel-optical-depth (POD) analysis and com- 
pared with forward-modeled simulations of the IGM to 
infer [C/H] versus redshift. 

The POD method yields a median [C/H]= -3.47 at a 
pathlength- weighted mean redshift of 3.1. This lies be- 
tween our estimate of [C/H]=-3.1 at z ~ 2.4 and [C/H]=- 
3.6 at z ~ 4.3. A simple linear interpolation between our 
two redshift points would yield [C/H]=-3.3 at z = 3.1, 
about 50% higher than the POD estimate. Schaye et 
al. find no statistically significant evidence of redshift 
evolution, in contrast to our result. 

Some of this discrepancy may reflect a real difference in 
the measurements and ionization corrections. But there 
are several factors to consider when comparing the two. 
First, while we and Schaye have both employed the up- 
dated Haardt & Madau model of the UV background 
with galaxies included, we have softened our spectrum 
above 10 Ryd at z = 4.3, to avoid a rise in the X-ray 
background that would be inconsistent with the observed 
decline in the space density of AGN. Schaye et al noted 
that in experimenting with models softened above 4 Ryd 
(simulating He II reionization), a transition in the middle 
of their redshift range would result in a positive detection 
of evolution . For reference, when we use th e updated cal- 
culations of iFaucher-Giguere et al.l (|2009T ) which are not 
artificially tuned but which incorporate He II ionization 
into the mean calculation, we measure an even stronger 
evolution by about 0.3 dex. 

A second point is that the present study (by design) 
has its pathlength weighted at the two endpoints of the 
redshift interval where we have detected a time deriva- 
tive in the metallicity. The Schaye et al sample has its 
pathlength weighted toward the middle of its redshift 
range (around z ~ 3.0) to maximize signal-to-noise, and 
therefore has somewhat less leverage for measuring evo- 
lutionary trends. For example, of the 19 QSOs in their 



sample, only 2 have pathlength above z = 3.6, and only 
one has pathlength above z = 4.0. Indeed, the relative 
lack of high redshift, high resolution data in the Schaye 
sample was a motivating factor for obtaining our MIKE 
data set. 

This effect manifests partly as a reduced signal-to-noise 
for measuring evolution. However, it is also the case that 
the lowest densities probed in any large sample spanning 
Az > 1 will be found at higher redshift, because of the 
reduced ionization of the IGM. Likewise the highest red- 
shifts are best for measuring low [C/H] because of the 
large fc iv at earlier times. Since density, abundance, 
and redshift can be artificially correlated, a large regres- 
sion analysis like the POD technique may encounter de- 
generacies along these basis vectors. The POD analysis 
projects a strong signal along the density axis with weak 
redshift evolution. By studying slices of samples with in- 
dividual absorbers, we favor a slightly weaker correlation 
with density, and a stronger evolution with redshift. 

Finally, we note that simulations of metal enrich- 
ment in the IGM ha ve improved substantia l ly ove r the 
last several years. lOppenheimer fc Pavel (|2006l ) and 
iCen fc Chisaril ([20101 ) have performed extensive stud- 
ies of metal enrichment in cosmological simulations 
with particular attention to carbon and the C iv ion. 
These represent imp roved updates to the early work of 
lAguirre et ail (|2001l ). They find that the total budget of 
heavy elements in the IGM increases by a factor of 2.7, or 
0.4 dex over this range. This is very closely in line with 
the change that we observe, indicating that wind mod- 
els motivated by local observations-but implemented at 
high rcdshift-provide a good match to the trends seen in 
the carbon abundance. 

8. CONCLUSIONS 

We have derived the distribution function of carbon 
abundance in Ly-a forest clouds at z ~ 4.3, using a set 
high signal-to-noise ratio spectra taken with the MIKE 
echelle spectrograph on the Magellan Clay telescope. 
The C iv column density or its upper limit was mea- 
sured for an H I -selected sample of 131 discrete absorbers 
with Ah j > 10 14 5 cm -2 , corresponding to p/p > 1.6. 
These measurements were converted into [C/H] abun- 
dances via density-dependent ionization corrections. The 
[C/H] distribution was then determined via the Kaplan- 
Meier product limit estimator for censored data. Our 
main findings are that: 

1. Over the range we can probe, the C iv distribution 
at z ~ 4.3 appears to be crudely lognormal, with a 
median of [C/H]=-3.55 and scatter of 0.8 dex. 

2. The median abundance at z ~ 4.3 is about 0.3-0.5 
dex lower than at z ~ 2.4, at fixed cosmic overden- 
sity p/p). The range quoted reflects differing as- 
sumptions about the UV ionizing background, and 
the gradient of abundance with density. 

3. We examined several sources of uncertainty in the 
measurements, including spatial variations in the 
radiation field and misestimates of the continuum. 
Variations in the background may contribute to 
scatter in abundance estimates at the ~ 0.2 dex 
level, but do not contribute a significant system- 
atic error. The same is true for continuum errors. 



19 



If systematic errors are present they would tend 
to enhance, rather than diminish, the evolution- 
ary signal. The one exception to this rule is the 
proposed "sawtooth" shaped attenuation of the 3-4 
Ryd background from He II absorption. Absorption 
at the levels indicated in iMadau fc Haardtl (2009) 
lead to a systematic increase of 0.06 dex in our 
abundances. A substantially stronger absorption 
signature could weaken the detected evolution sig- 
nal, but this may also cause conflict with measure- 
ments of the Si iv /C iv ratio in the limited cases 
where this could be measured. 

4. The total carbon contribution to closure density in 
our fiducial model is fie ~ 2.7 x 10~ 8 . The same 
quantity calculated at z ~ 2 is 1.7 times larger, 
implying that roughly half of the heavy elements 
seen in the Ly-a forest at z ~ 2.4 were distributed 
into the IGM in the 1.3 Gyr of the z ~ 3—4 interval. 

5. The mass flux of carbon into the IGM needed to 
sustain this growth is comparable in magnitude to 
the carbon yield inferred from newly formed stellar 
populations at these redshifts. We estimate that on 
global scales (and within very large uncertainties) , 
the metal feedback rate from galaxies is around 
~ 30% of the star formation rate. This large rate is 
consistent with observations of star forming galax- 
ies in the local universe. 

The main result of the paper — our detection of a de- 
cline in abundance at higher redshift — shows that at least 
some of the enrichment of the IGM is taking place at the 
time when we observe it. High redshift C iv observa- 
tions show that some metals already existed at z ~ 5.5, 
but a substantial fraction of the heavy elements seen at 
z 2.4 must have been ejected from galaxies we can 
readily observe. 

The transport of metals from galaxies to the IGM is 
one of the strongest lines of evidence supporting the no- 
tion that feedback is a crucial element of galaxy forma- 
tion. This has long been inferred indirectly, based on 
observations and models of galaxy properties in the local 
universe. Our measurements suggest that observations at 
high redshift may start to place meaningful constraints 
on the feedback process during the epoch when it was 
actually taking place. 
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<1.521730e+12 






4.427306 





000053 


23 


08 


4.50 


14 


588 


0.182 


3.175900c+12 


5 


549000c+ll 


4.387300 





000023 


27 


53 


0.86 


15 


016 


0.030 


1.453000e+12 


3 


534400c+ll 


4.365353 





000028 


26 


82 


2.47 


14 


979 


0.056 


4.482700c+12 


1 


404300e+12 


4.363416 





000027 


30 


34 


2.41 


15 


810 


0.017 


6.620800e+12 


1 


705800e+12 


4.362077 





000001 


41 


84 


2.39 


15 


129 


0.030 


6.537800c+12 


1 


692900e+12 


4.358888 





000098 


49 


68 


2.09 


16 


601 


0.086 


1.406300c+13 


1 


466300e+12 


4.357327 





000045 


37 


50 


4.37 


15 


950 


0.039 


1.499900c+13 


2 


223700e+12 


4.355646 





000039 


42 


65 


1.00 


15 


898 


0.038 


9.965700c+12 


9 


962500c+ll 


4.347564 





000036 


22 


04 


2.38 


15 


015 


0.057 


9.183500c+12 


1 


373000e+12 


4.349267 





000027 


39 


98 


4.95 


16 


390 


0.219 


9.769600c+12 


1 


266500c+12 


4.345974 





000007 


25 


19 


0.39 


14 


877 


0.025 


<2.154060e+12 






4.321407 





000127 


49 


20 


3.67 


15 


572 


0.083 


1.916900c+13 


1 


270200c+12 


4.320024 





000000 


43 


28 


2.31 


16 


145 


0.091 


6.301600c+13 


2 


671100c+12 


4.317857 





000028 


23 


37 


3.56 


14 


758 


0.176 


<1.935430e+12 






4.312235 





000006 


24 


84 


1.16 


14 


610 


0.073 


<1.957310e+12 






4.283492 





000007 


22 


66 


0.66 


15 


313 


0.056 


<9.478000c+ll 






4.279064 





000460 


29 


93 


4.83 


14 


599 


1.884 


<7.070800c+ll 






4.275294 





000028 


29 


63 


1.73 


15 


189 


0.053 


<9.456400c+ll 






4.270005 





000004 


26 


80 


0.52 


15 


552 


0.057 


<1.662370e+12 






4.266536 





000093 


20 


03 


7.51 


14 


912 


0.365 


<1.636960e+12 






4.265502 





000250 


17 


01 


12.58 


14 


652 


0.484 


<2.184220e+12 






4.262394 





000010 


25 


89 


1.45 


14 


637 


0.045 


<5.307900e+12 






4.251741 





000005 


18 


25 


0.83 


14 


509 


0.067 


<1.670080e+12 






4.248899 





000071 


38 


35 


2.15 


15 


081 


0.057 


<1.354510c+12 






4.246332 





000010 


25 


62 


0.54 


14 


821 


0.017 


<2.489940c+12 






4.215746 





000026 


37 


14 


1.10 


15 


069 


0.023 


<1.363900e+12 






4.214093 





000054 


35 


74 


2.14 


14 


524 


0.066 


3.027200c+12 


4 


501900c+ll 


4.211229 





000044 


43 


14 


3.19 


15 


521 


0.095 


2.055800c+13 


2 


462300e+12 


4.207951 





000034 


38 


72 


2.55 


15 


050 


0.040 


1.854900c+12 


5 


446300c+ll 


4.203494 





000011 


20 


59 


0.91 


14 


538 


0.018 


<2.748370e+12 






4.195487 





000021 


38 


94 


3.37 


16 


121 


0.239 


1.298400c+13 


1 


459910c+ll 


4.187106 





000008 


60 


09 


1.08 


15 


146 


0.032 


<2.299030e+12 






4.181189 





000317 


43 


20 


9.22 


14 


605 


0.357 


<2.430070c+12 






4.179194 





000020 


31 


28 


1.32 


15 


106 


0.025 


<3.316110e+12 






4.159921 





000068 




DZ 


6.29 


14 


940 


0.046 


<6.238000e+ll 






4.157869 





000161 


57 


96 


19.39 


14 


547 


0.127 


<7.426000e+ll 






4.15DZ32 





000063 


22 


15 


5.40 


15 


194 


0.196 


2.136000e+12 


3 


455100e+ll 


4.150681 





000009 


24 


02 


0.83 


15 


046 


0.070 


6.060600c+ll 


1 


882100c+ll 


4.128559 





000006 


25 


43 


0.87 


15 


121 


0.089 


<9.137800c+ll 






4.084159 





000007 


28 


10 


1.87 


15 


369 


0.155 


6.524400c+12 


1 


244600c+12 


4.065986 





000011 


20 


46 


0.76 


14 


752 


0.033 


2.401500c+12 


7 


137400c+ll 


4.061522 





000054 


20 


45 


2.01 


14 


715 


0.109 


<1.363140e+12 






4.060768 





000109 


29 


66 


13.62 


14 


517 


0.202 


<7.458400c+ll 






4.058138 





000049 


41 


92 


6.98 


14 


858 


0.053 


<2.221440e+12 






4.056497 





000082 


35 


69 


11.98 


14 


786 


0.250 


2.833900c+12 


5 


575700c+ll 


4.048319 





000090 


45 


59 


9.39 


16 


136 


0.377 


1.334600c+13 


1 


830900c+12 


4.039261 





000006 


29 


48 


0.66 


14 


956 


0.034 


<3.718200e+12 






4.024303 





000001 


20 


00 


1.00 


14 


710 


0.168 


<1.545500e+12 







Table 2 

Line Sample H I and C IV Measurements (Continued) 





ZH I 




dz 


bn i 


db 




dN H1 


Wciv 


dNciv 
















BR0418-5 


726 






4. 


353709 





000011 


29 


.37 


0.63 


15.409 


0.029 


<4.045600c+12 




4. 


352097 





.000025 


17 


.44 


1.29 


14.720 


0.035 


<3.988130e+12 




4. 


343351 





.000021 


15 


.91 


1.56 


15.818 


0.162 


4.201200c+12 


6.709000c+ll 


4. 


326431 





.000022 


42. 


.55 


1.10 


14.991 


0.023 


<3.344200e+12 




4. 


320692 





.000033 


27 


.35 


1.15 


14.608 


0.038 


1.187000c+12 


3.475000e+ll 


4. 


296004 





.000021 


23 


.62 


1.35 


16.636 


0.207 


<2.976100e+12 




4. 


294459 





.000045 


28 


.64 


2.03 


15.140 


0.046 


<2.487100e+12 




4. 


279348 





.000007 


13 


.89 


1.78 


14.681 


0.293 


<1.675300e+12 




4. 


274486 





.000062 


28 


.09 


1.65 


15.243 


0.102 


<1.150750e+12 




4. 


271221 





.000130 


46 


.95 


5.46 


14.528 


0.083 


<2.717920e+12 




4. 


259756 





.000013 


19 


.91 


0.62 


14.636 


0.025 


<1.717482e+12 




4. 


249553 





.000010 


28 


.42 


1.36 


14.812 


0.020 


<4.878660c+12 




4. 


246544 





.000013 


39 


.49 


1.43 


15.539 


0.064 


<1.117540e+12 




1. 


236839 





.000185 


27 


.35 


6.57 


16.190 


0.421 


1.896900c+13 


5.357000e+12 


4. 


238187 





.001232 


58 


.60 


130.50 


15.636 


1.060 


<3.061680e+13 




4. 


239286 





.000301 


29 


.97 


23.48 


15.664 


0.703 


2.149500c+13 


2.059100c+12 


1. 


221174 





.000058 


43 


.81 


5.60 


15.847 


0.081 


<1.110400e+14 




4. 


219525 





.000039 


15 


.52 


8.52 


17.091 


2.256 


4.352500c+13 


2.246000e+12 


4. 


218226 





.000215 


45 


.33 


16.66 


15.324 


0.315 


<2.476850e+12 




4. 


215030 





.000048 


42 


.13 


3.38 


14.688 


0.044 


<2.010750e+12 




4. 


209963 





.000034 


29 


.78 


1.16 


15.354 


0.050 


1.933300c+12 


4.152700c+ll 


4. 


204677 





.000063 


21 


.10 


3.52 


15.364 


0.632 


<7.646100e+12 




4. 


199868 





.000041 


24 


.82 


3.46 


15.269 


0.064 


1.329100c+13 


1.268400e+12 


4. 


198389 





.000330 


46 


.81 


20.02 


15.428 


0.257 


<3.428510e+12 




4. 


196863 





.000250 


50. 


.00 


7.49 


15.700 


0.164 


2.229900c+12 


7.301700c+ll 


4. 


176287 





.000007 


21 


.72 


1.53 


14.939 


0.154 


<2.710600c+12 




4. 


174260 





.000057 


36 


.30 


2.14 


14.627 


0.094 


<1.355270e+12 




4. 


159184 





.000009 


26 


.70 


1.64 


15.510 


0.192 


<2.615960e+12 




4. 


139674 





.000069 


38 


.80 


1.99 


15.483 


0.085 


<5.521110e+12 




4. 


130618 





.000315 


19 


.02 


5.50 


14.678 


0.898 


<1.310500e+12 




4. 


124536 





.000013 


34 


.26 


1.95 


15.057 


0.109 


<3.222700e+12 




4. 


120722 





.000009 


47 


.18 


1.94 


14.813 


0.044 


<3.442600e+12 




4. 


108080 





.000030 


59 


.86 


2.83 


15.315 


0.053 


<3.456100e+12 




4. 


094620 





.000008 


24 


.13 


0.78 


15.459 


0.086 


<2.806240e+12 




1. 


084853 





.000036 


29 


.98 


2.82 


14.871 


0.096 


<6.064000c+ll 




4. 


081859 





.000013 




no 

yy 


1.70 


15.108 


0.070 


<3.831100e+12 




4. 


079143 





.000011 


34 


.27 


2.57 


15.588 


0.113 


2.185100c+13 


1.434600c+12 


4. 


065194 





.000018 


30 


.88 


1.25 


15.028 


0.047 


<1.995340e+12 




4. 


063509 





.000061 


28 


.00 


2.88 


14.678 


0.060 


<2.679990e+12 




4. 


056032 





.000014 


33 


.74 


1.00 


15.077 


0.034 


3.049600c+12 


6.418300c+ll 


4. 


054350 





.000084 


27 


.37 


3.88 


14.561 


0.104 


<2.664230e+12 




4. 


053611 





.000115 


16 


.06 


11.71 


14.565 


0.231 


<3.405100e+12 




4. 


052538 





.000064 


25 


.00 


1.41 


16.167 


0.304 


<1.471000e+12 




4. 


019813 





.000011 


27 


.11 


1.56 


14.829 


0.099 


<2.463100e+12 




4. 


017804 





.000012 


35 


.88 


3.03 


15.188 


0.156 


<3.272670e+12 




4. 


013256 





.000022 


40 


.30 


4.53 


15.098 


0.159 


9.325000c+12 


1.006000c+12 


4. 


008585 





.000011 


17 


.36 


2.86 


14.616 


0.321 


<1.370500e+12 




1. 


006144 





.000013 


27 


.81 


2.45 


14.742 


0.134 


<5.012040e+12 




4. 


003678 





.000011 


30 


.28 


2.37 


15.187 


0.184 


<3.619320c+12 





Table 2 

Line Sample H I and C IV Measurements (Continued) 



ZH I 




dz 


I 


db 




dN H1 




dN c iv 


BR0714-6449 


4.397271 





.000008 


28 


.22 


0.71 


15.140 


0.054 


<2.065630c+12 




4.394346 





,000026 


46 


.01 


1.19 


15.321 


0.026 


6.065800c+12 


1.274700c+12 


4.391425 





.000034 


12 


.99 


5.94 


16.467 


1.750 


<1.431670c+12 




4.388564 





.003111 


47 


10 


64.55 


15.540 


3.027 


1.953600c+13 


1.720900c+12 


4.383436 





.000032 


45 


08 


2.72 


15.412 


0.068 


<1.610860c+12 




4.381303 





.000027 


30 


.27 


1.45 


15.343 


0.068 


<1.665610c+12 




4.376600 





.000032 


48 


01 


1.89 


14.808 


0.028 


<2.499250c+12 




4.371056 





.000055 


25 


.82 


2.18 


15.598 


0.196 


<1.634500c+12 




4.363155 





.000024 


48 


.49 


1.19 


15.818 


0.052 


<5.346400e+12 




4.359026 





.000014 


38 


34 


1.97 


14.632 


0.067 


<3.695800e+12 




4.345805 





.000011 


28 


01 


1.79 


14.602 


0.098 


3.804600c+12 


1.078000c+12 


4.342708 





.000038 


29 


.64 


3.74 


14.516 


0.101 


<2.709610e+12 




4.339696 





.000008 


16 


.83 


0.44 


14.933 


0.057 


<2.503270c+12 




4.333522 





.000012 


12 


46 


2.88 


15.394 


0.873 


<3.287530c+12 




4.330081 





.000011 


32 


.24 


0.89 


15.838 


0.083 


<2.529160e+12 




4.326107 





.000015 


30 


.21 


1.60 


15.421 


0.087 


<2.999980e+12 




4.323170 





.000023 


32 


.99 


1.69 


15.390 


0.080 


<3.884960e+12 




4.318676 





.000013 


21 


.08 


0.75 


14.726 


0.035 


<7.277800e+12 




4.290688 





.000008 


18 


.97 


1.76 


15.522 


0.344 


<3.648910e+12 




4.272618 





.000009 


19 


92 


1.06 


15.098 


0.133 


2.932500c+12 


7.360300c+ll 


4.269980 





.000040 


21 


,24 


1.47 


14.926 


0.100 


<8.070100c+ll 




4.261784 





.000452 


38 


76 


24.13 


14.603 


0.358 


<4.188790c+12 




4.243142 





.000009 


18 


,81 


1.23 


15.648 


0.190 


<2.221580e+12 




4.233165 





.000021 


20 


.75 


1.08 


15.199 


0.091 


<3.497500e+12 




4.213589 





.000048 


32. 


81 


3.06 


15.227 


0.049 


<2.230290e+12 




4.200533 





.000030 


44. 


15 


2.22 


14.712 


0.032 


<2.774050e+12 




4.197665 





.000312 


42 


,29 


6.58 


16.316 


0.700 


1.792100c+14 


2.178100c+13 


4.187770 





.000034 


26 


30 


1.45 


14.791 


0.047 


<2.836570e+12 




4.185000 





.000156 


70 


00 


13.32 


15.023 


0.073 


<9.562500e+12 




4.183234 





.000033 


33 


,25 


2.95 


15.488 


0.065 


7.748100c+12 


8.997800c+ll 


4.181524 





.000024 


23 


,86 


2.08 


15.039 


0.112 


<1.065910e+12 




4.169388 





.000035 


41 


.54 


1.49 


15.055 


0.034 


<2.024560e+12 




4.159509 





.000039 


22 


.68 


4.96 


16.186 


1.037 


<5.225800c+12 




4.155559 





.000783 


25 


.00 


23.55 


15.071 


1.210 


<2.646090c+12 




4.150633 





.000462 


31 


29 


18.51 


15.541 


1.000 


5.041800c+12 


1.463100c+12 


4.149312 





.000054 


20 


76 


3.25 


15.040 


0.125 


<1.921510c+12 




4.143332 





.000050 


48 


19 


3.01 


15.036 


0.038 


2.864000c+12 


6.603200c+ll 


4.139893 





.000042 


70. 


.00 


5.18 


16.185 


0.172 


<1.117390e+12 




4.128364 





.000017 


30 


.32 


2.31 


15.949 


0.236 


<1.328410c+12 




4.122314 





.000019 


18 


16 


3.56 


16.404 


0.953 


<2.115940c+12 





